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p^. Abstract 

I We study the possibility of producing a new kind of nuclear systems which in addition to 

Os . _ _ 

ordinary nucleons contain a few antibaryons {B = p, A, etc.). The properties of such systems 

are described within the relativistic mean-field model by employing G-parity transformed inter- 



> 

! actions for antibaryons. Calculations are first done for infinite systems and then for finite nuclei 

' from ^He to ^'^^Pb. It is demonstrated that the presence of a real antibaryon leads to a strong 

■ rearrangement of a target nucleus resulting in a significant increase of its binding energy and 
local compression. Noticeable effects remain even after the antibaryon coupling constants are 

Q I reduced by factor 3 — 4 compared to G-parity motivated values. We have performed detailed 

^ ' calculations of the antibaryon annihilation rates in the nuclear environment by applying a kinetic 



X 



approach. It is shown that due to significant reduction of the reaction Q-values, the in-medium 
annihilation rates should be strongly suppressed leading to relatively long-lived antibaryon- 
nucleus systems. Multi-nucleon annihilation channels are analyzed too. We have also estimated 
formation probabilities of bound B + A systems in pA reactions and have found that their ob- 
servation will be feasible at the future GSI antiproton facility. Several observable signatures are 
proposed. The possibility of producing multi-quark-antiquark clusters is discussed. 

PACS numbers: 25.43.+t, 21.10.-k, 21.30.Fe, 21.80.+a 
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I. INTRODUCTION 



Antibaryons are extremely interesting particles for nuclear physics. They are building 
blocks of antimatter which can be produced in the laboratory. When in close contact 
baryons and antibaryons promptly annihilate each other producing mesons in the final 
state. These reactions have been actively studied in 80's and 90's, most extensively at 
LEAR (see reviews [1, 2]). In particular, many attempts have been made to find A^A^ 
bound states with mass close to the threshold [3], but up to now no clear evidence was 
found (see, however, the recent paper [4]). 

In contrast to elementary BB interactions, much less is known about antibaryon- 
nuclear interactions. Main information in this case is provided by antiprotonic atoms 
and scattering data. However, due to strong absorption this information is limited 
to a far periphery of the nuclear density distribution. As follows from the analysis of 
Refs. [5, 6], the real part of the antiproton optical potential in nuclei might be as large 
as —(200 300) MeV with uncertainty reaching 100% in the deep interior. The imaginary 
part is also quite uncertain in this region. 

On the other hand, many interesting predictions concerning the antibaryon behavior 
in nuclear medium have been made. In particular, appearance of antinucleon bound 
states in nuclei is one of the most popular conjectures. This possibility was first studied 
in early 80's [5, 7-9] using antinucleon optical potentials consistent with the p-atomic 
data. However, much deeper antinucleon potentials have been obtained [10, 11] on the 
basis of relativistic nuclear models [12, 13]. They predict a large number of deeply bound 
antinucleon states in nuclei [14, 15]. Several observable signatures of such states have been 
discussed in Refs. [14, 16-18]. Besides, it was demonstrated in Ref. [19] that in-medium 
reduction of the antibaryon masses may lead to enhanced yields of baryon-antibaryon 
pairs in relativistic heavy- ion collisions. This mechanism was also used in Refs. [20, 21] 
to explain subthreshold antibaryon production in pA and AA collisions. 

The Relativistic Mean-Field (RMF) models [22, 23] are widely used now for describing 
nuclear matter and finite nuclei. Within this approach nucleons are described by the Dirac 
equation coupled to scalar and vector meson fields. Scalar S and vector V potentials 
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generated by these fields modify the spectrum of the Dirac equation in homogeneous 
isospin-symmetric nuclear matter as follows 

E^{p)^V±^/(^^^i:^^'^W+^ : (1) 

where mjv and p are the vacuum mass and 3-momentum of a nucleon, respectively. The 
+ sign corresponds to nucleons with positive energy 

E^{p) = E+{p) + V(m;v--5)2 + p2 , (2) 

and the — sign corresponds to antinucleons with energy 

^^(p) = -^-(-p) = -V+ ^(m^-Sy + p^. (3) 

Changing sign of the vector potential for antinucleons is exactly what is expected from the 
G-parity transformation of the nucleon potential. As follows from Eq. (1), the spectrum 
of single-particle states of the Dirac equation in nuclear environment is modified in two 
ways. First, the mass gap between positive- and negative-energy states, 2(mjv — -S"), is 
reduced due to the scalar potential and, second, all states are shifted upwards due to the 
vector potential. 

It is well known from nuclear phenomenology that a good description of nuclear ground 
state is achieved with S ~ 350 MeV and V ~ 300 MeV so that the net potential for 
slow nucleons is V — S* ~ —50 MeV. Using the same values one obtains for antinucleons 
a very deep potential, — — 5* ~ —650 MeV. Such a potential would produce many 
strongly bound states in the Dirac sea. However, when these states are occupied they are 
hidden from the direct observation. Only creating a hole in this sea, i.e. inserting a real 
antibaryon into the nucleus, produces an observable effect. If this picture is correct one 
can expect the existence of strongly bound states of antinucleons in nuclei. Qualitatively 
the same conclusions can be made about antihyperons (A, E, 5, . . .) which are coupled to 
the meson fields generated by nucleons. 

Of course, these bound antibaryon-nuclear systems will have limited life times because 
of unavoidable annihilation of antibaryons. The naive estimates based on the vacuum 
annihilation cross sections would predict life times of the order of 1 fm/c. However, due 
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to strong reduction of phase space available for annihilation in medium, these life times 
can be much longer. When calculating the structure of bound antibaryon-nuclear systems 
we first consider them as stationary objects. Afterwards their life times are estimated on 
the basis of kinetic approach taking into account in-medium effects. 

In our previous paper [24] we have performed self-consistent calculations of bound 
antibaryon-nuclear systems which take into account rearrangement of nuclear structure 
due to the presence of a real antibaryon. We have found not only a significant increase in 
the binding energy but also a strong local compression of nuclei induced by the antibaryon. 
The calculations were mainly done for antiprotons with potentials obtained by G-parity 
transformation. In the present work we extend our study in several directions. We 
report on new results for antiprotons and antilambdas for wider range of target nuclei. 
Calculations for reduced couplings of antibaryons to the meson fields are done too in order 
to simulate effects which are missing in a simple mean-field approximation. We have also 
performed a detailed study of the annihilation processes which determine the life time 
of the antibaryon-nuclear systems. The possibility of a delayed annihilation due to the 
reduction of the available phase space is demonstrated by direct calculations. Different 
scenarios of producing bound antibaryon-nuclear systems by using antiproton beams are 
considered and their observable signatures are discussed. In particular, we point out the 
possibility of inertial compression of nuclei and deconfinement phase transition induced 
by antibaryons. 

Our paper is organized as follows. In Sect II we describe the formalism used in calcula- 
tions of antibaryon-nuclear systems and briefly explain the numerical procedure applied 
for solving the RMF equations. The limiting case of infinite matter with admixture of 
antibaryons is studied in Sect. III. The numerical results concerning bound antibaryon- 
nuclear systems are presented in Sect. IV. The problem of antibaryon annihilation both in 
infinite and finite nuclear systems is discussed in Sect. V. Probabilities of creating nuclei 
with bound antiprotons and antilambdas in pA collisions are estimated in Sect. VI. Pos- 
sible observable signatures of such systems are considered in Sect. VII. Finally, in Sect. 
VIII we summarize our results, formulate several open questions and tasks for future. 
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II. THEORETICAL FRAMEWORK 

To study antibaryon-nuclear bound states we use different versions of the relativistic 
mean-field (RMF) model. This approach has several advantages: a) its effective de- 
grees of freedom are adequate for the description of nuclear many-body systems, b) it 
describes nuclear ground-state observables with high accuracy, c) antinucleons are nat- 
urally included in the theory through the nucleonic Dirac equation, d) its generalization 
for antihyperons is straightforward . 

There remain, however, uncertainties in its application to antibaryons. We list only 
some of them: a) the G-parity symmetry is not necessarily valid in a many-body system 
(we shall discuss this issue later), b) the coupling constants of the model as well as its 
functional form are determined for nuclear ground states and extrapolation to higher 
densities is somewhat uncertain, c) it is believed that the RMF model can be considered 
as an effective field theory only at low enough energies, A£' < 1 GeV. On the other hand, 
inserting an antibaryon into a nucleus delivers an excitation energy of about 2 GeV, which 
might be beyond the applicability domain of such models. 

Despite of these uncertainties, the RMF approach allows an exploratory study of new 
types of finite nuclear systems. The RMF model (for reviews, see [22, 23]) is formulated 
in terms of a covariant Lagrangian density of nucleons and mesons. We modify the 
usual approach by introducing real antibaryons in addition to protons and neutrons. The 
Lagrangian density L reads i^ — c— 1) 



j=N,B 



+ 




— -F F^" 



(4) 
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Here the degrees of freedom are nucleons (iV), antinucleons (B = N) or antihyper- 
ons {B = A, S,...), the isoscalar a {J'" = 0^) and uj {J'^ = 1^) mesons, the isovector 
p ( J'^ = 1~) meson as well as the Coulomb field A^. The field tensors G^j^ of the vector 
fields {G — cu, p,A) are defined as 

G/xu = d^Gv — dvG^ . (5) 

Arrows refer to the isospin space. The quantities rrij in Eq. (4) denote the vacuum masses 
of nucleons and antibaryons, m^^ and rUp are, respectively, the vacuum masses of uj and p 
mesons.The a meson mass and the coupling constants gij are adjustable parameters 
which are found by fitting the observed data on medium and heavy nuclei (see below) . 

Our approach employs two common approximations, the mean-field and the no-sea 
approximations. The latter implies that we consider explicitly only the valence Fermi and 
Dirac sea states. On the other hand, the filled Dirac sea states are included only imphcitly 
via nonlinear terms of scalar and vector potentials. According to recent work [25], the 
vacuum contribution and the related counter terms have a structure similar to the meson 
field terms. Implementing a real antinucleon inside a target nucleus means appearance of 
a hole in the nucleonic Dirac sea. The mean-field approximation consists in replacing the 
meson field operators and the corresponding source currents by their expectation values. 
This leads to classical mean fields, for example: 

a^{a), (6) 

Below we treat the RMF model in Hartree approximation i.e. neglect all exchange terms 
arising from antisymmetrization of single-particle states. 

Equation of motion for all meson and fermion fields are obtained from the Euler- 

Lagrange equations 

d f dC \ dC ^ _ ^ ^ 

U(%/^^^) J ' % ^ ' ,.^N,B,a,.,pA. (7) 

In the equilibrium (static) case this gives the stationary Dirac equations for (anti)baryons 

e"^^"^ = [-i otV + P {mj -Sj) + Vj\ i^f . (8) 
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Here j — N,B and a denotes various single-particle valence states with the wave func- 
tions ip'^ and energies . The scalar and vector potentials acting on (anti)baryons are 
defined as 



Sj = g^jcr, (9) 
V, = g^ju;' + g,,plT, + e,A'^-^ . (10) 



Since we consider static systems with even numbers of neutrons and protons, only 
time-like components of vector fields give nonzero contributions in the mean-field ap- 
proximation. Disregarding the mixing of proton and neutron states we retain only the pa 
components of the p-meson field. We also assume that the time-reversal invariance is 
valid even in the presence of an unpaired particle like an antiproton. 

Equations of motion for the meson fields read 

{-/^ + ml + ha + c(t'^)(t = ^g^jPsj, (11) 

j 

{-A + ml + du;o)u;o = ^QwjPj, (12) 

j 

(-A + m2)p0 = Y.9P3PI3^ (13) 

3 

AAo = Y.^jPQj (14) 

3 

The scalar, vector, isovector and charge densities are defined as 

Psj = (V'jV'j), (15) 
P3 = (V'M> , (16) 

Pi3 = {i']ni^3) > (17) 

PQ3 = ^(^^1(1 + ^3)^^,), (18) 

where angular brackets mean averaging over the ground state wave function. 

To investigate sensitivity of the results to the model parameters, we have consid- 
ered several versions of the RMF model, namely the NL3 [26], NLZ [27], NLZ2 [28] and 
TMl [29] models. The corresponding parameter sets are fisted in Table I. They were 



7 



TABLE I: The paramotor sots of the RMF models used in this paper. 





NLZ 


NLZ2 


NL3 


TMl 


mjv (MeV) 


938.9 


938.9 


939.0 


938.0 


nia (MeV) 


488.67 


493.150 


508.194 


511.198 


(MeV) 


780.0 


780.0 


782.501 


783.0 


nip (MeV) 


763.0 


763.0 


763.0 


763.0 


9aN 


10.0553 


10.1369 


10.217 


10.0289 


guiN 


12.9086 


12.9084 


12.868 


12.6139 


9pN 


4.84944 


4.55627 


4.4740 


4.6322 


b (fm-^) 


13.5072 


13.7561 


10.431 


7.2325 


c 


-40.2243 


-41.4013 


-28.885 


0.6183 


d 











71.305 



found by fitting binding energies and radii of spherical nuclei from ^^O (not included in 
the TMl fit) to Pb isotopes [67]. 

Regarding the antibaryon couplings, there is no reliable information suitable for high 
density nuclear matter. In this situation, as a starting point, we choose antibaryon-meson 
coupling constants motivated by the G-parity transformation. It is analogous to the 
ordinary parity transformation in the configurational space, which inverts the directions 
of 3-vectors. The G-parity transformation is defined as the combination of the charge 
conjugation and the 180° rotation around the second axis of the isospin space [30]. As 
known, a and p mesons and the Goulomb field have positive G-parity, while cu meson 
has negative G-parity. Therefore, applying the G-parity transformation to the baryon 
potentials (9)-(10) one obtains the corresponding potentials for antibaryons. The results 
of this transformation can be formally expressed by the following relations between the 
coupling constants: 

9aN = 9aN , 9,^17 = -9loN , 9 pN = 9pN ■ (19) 

If these relations are valid, one can make two conclusions. First, equal scalar potentials 
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lead to equal effective masses for nucleons and antinucleons. Second, the vector poten- 
tials have opposite signs for these particles. This is in agreement with the conclusion 
made by considering positive- and negative-energy solutions of the Dirac equation (see 
Eqs. (2)-(3)). 

The simple consideration based on the G-parity transformation of mean meson fields 
is certainly an idealization. There are several effects in many-body systems which can, 
in principle, significantly distort this picture. Here we mention only two of them. First, 
inclusion of exchange (Fock) terms leads generally to smaller magnitudes of scalar and vec- 
tor potentials for antinucleons in dense baryon-rich matter as compared to nucleons [31]. 
Second, the expressions (9)-(10) correspond to the tadpole (single bubble) diagrams as- 
sociated with the (7,0; and p meson exchanges. However, for antibaryons one should 
consider also the contribution of more complicated multi-meson diagrams originating 
from annihilation channels in the intermediate state. By using dispersion relations, the 
corresponding contribution to the real part of the antibaryon self energy can be expressed 
through the BN annihilation cross section. The calculations of Refs. [20, 21] show that 
such a contribution can reach 100 — 150 MeV for slow antinucleons at normal nuclear 
density. 

Taking into account all these uncertainties, in our calculations we consider not only 
the G-parity motivated couplings (19), but also the reduced antibaryon couplings 

9aN = ^9tN , g^N = -igujN , Q pjj = ^QpN (20) 

with the same reduction parameter ^ for interactions with cr, uj and p meson fields. 
By choosing ^ in the interval < ^ < 1 we can investigate all possibilities from max- 
imally strong antibaryon couplings to noninter acting antibaryons. In calculations of 
antihyperon-nuclear system we use the coupling constants motivated by the SU(3) flavor 

symmetry, namely assume that g^j^ = ^g^j^, g^-^ = ^g^-^ . 

By using Eq. (20) we can rewrite the source terms in Eqs. (11)-(12) in a more trans- 



9 



parent form 

^ QajPSj = QaN {PSN + ^Psn) : (21) 

j 

Yl dojjPj = 9ojN (pn - CPn) ■ (22) 
j 

One can infer from these relations that presence of antibaryons in nuclear matter leads to 
increase of the scalar (attractive) potential and decrease of the vector (repulsive) poten- 
tial acting on surrounding nucleons. Unlike some previous works, we take into account 

the rearrangement of nuclear structure due to the presence of a real antibaryon This 
leads to enhanced binding and additional compression of nuclei as was first demonstrated 
in Ref. [24]. 

A few remarks about the numerical procedures are in place here. Calculations are per- 
formed for nuclear systems obeying the axial and reflection symmetries. We replace the 
Dirac equations (8) for single-particle wave functions V'" by the effective Schrodinger equa- 
tions for their upper components. Introducing explicitly the upper {ip+) and lower (</'-) 
spinor components of the wave function and omitting indices j and a , we rewrite Eqs. (8) 
in the form 

i}+\_lm-S+V (TP W 

ap -m + S + V y I "0- 
Here <t denote the spin Pauli matrices and p — — iV . Prom the lower part of Eq. (23) 
one gets 

11)^ = {e + m- S -V)-'^{ap)%lj+. (24) 
Substituting this expression into the upper part leads to the following equation for ■^^ 

eV^+ = hxl^+ , (25) 

where h is the effective, energy-dependent Hamiltonian 

h = (ap) (e + m-S-V)-^ {ap) +m-S + V. (26) 

Eqs. (25) for nucleon and antibaryon wave functions are solved by iterations, using 
the damped gradient method suggested in Ref. [32]. At each iteration step we use the 
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(23) 



potentials S, V determined by solving Eqs. (11)-(14) for mean meson and Coulomb fields 
at the preceding step. Numerical calculations are carried out on a spatial grid with 
equidistant grid points. The Fourier transformation of fermion and meson fields is used 
to evaluate spatial derivatives as matrix multiplication in the momentum space. This 
allows us to use a rougher grid with the cell size ~ 0.7 fm, as compared to much finer 
grids required by the finite difference schemes. Most calculations below were carried out 
with the cell size of 0.3 fm which guarantees the the numerical accuracy in binding energies 
and density profiles better than 0.5%. For nucleons we implement the paring correlations 
using the BCS model with a 6 pairing force and a smooth cut-off given by the Fermi 
distribution in single-particle energies [28]. The center of mass corrections are taken into 
account in the same way as for normal nuclei (for details see Ref. [28]). 

III. INFINITE MATTER WITH ADMIXTURE OF ANTIBARYONS 

To study qualitative effects due to the presence of antibaryons, let us consider first ho- 
mogeneous isospin-symmctric nucleonic matter with admixture of antibaryons (BN mat- 
ter). As compared to the general formalism developed in the preceding section, here 
we neglect the isospin-asymmetry and Coulomb effects assuming vanishing p-meson and 
electromagnetic fields. At zero temperature, disregarding medium polarization (density 
rearrangement) effects, we calculate thermodynamic functions of BN matter at fixed, 
spatially homogeneous densities of nucleons pn and antibaryons p^. In this case one can 
omit the spatial derivatives in Eqs. (11)-(14) and solve the Dirac equations (8) in the 
plane wave representation. For example, assuming that the G-parity concept is valid on 
the mean-field level, one gets the energy spectra of N, N states given by Eqs. (2)-(3) 
where S = q^nct, V = Qu^nuJo- 

Within the mean field approximation the occupation numbers of nucleons (j = N) and 
antibaryons {j — B) at zero temperature have the form of Fermi distributions Q{pFj — p) 
where Q{x) = (1 + sqna;)/2 and ppj = {Q'n'^ Pj /i^jY^^ is the Fermi momentum of corre- 
sponding particles (here Vj denote their spin-isospin degeneracy factors). Applying the 
standard formalism [22] one can obtain the energy density e — E/V = directly from 
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the Lagrangian (4) with the result 

e = Yl (S'" + 9o.jPjU;o) + + ^a' + ^a'- \mlul - ju'^ . (27) 



j=N,B 

Here e^'" is the " kinetic" part of the energy density 

PFj 



I dpp\/mf+p\ (28) 







where m* is the effective mass of j-th particles, related to the scalar meson field as 

m* = ruj - g^ja . (29) 
After calculating the integral in Eq. (28) we get the formula 

where ^(2;) is a dimensionless function 

1 2 4 

*(z) = 4 J dtt^VWTz^^VlT^(^l + ^^ -J^^ 


Equations for meson fields can be obtained by minimizing the energy density with 
respect to a and loq . Prom Eqs. (27)-(29) one gets 



^^±i. (31) 



{ml + ba + ca^)a ^ ^ g^jpsj : (32) 

j=N,B 

{ml + douDuo^ gu,jPj, (33) 

j=N,B 

where psj = de/drrij coincides with the scalar density defined in preceding section. The 
explicit expression for psj has the form 

Psj = iTi / dpp' ; ' = M , 34 







where 



*(,),_Li = vr+i^ (35) 
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FIG. 1: Energy per particle of homogeneous NN matter, e = E/{Nn + Nj^) , as a function of 
nucleon density (in units of pQ = 0.15 fm"^) at fixed % = pj^/pN- Upper (lower) panel shows 
results of the TMl (NLZ) calculation. Points A correspond to minima of e at x = 1/16 . 
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Using (33) one can rewrite Eq. (27) in a simpler form 

e = + + + + ^^''^^ + ■ ^^^^ 

3=N,B 

Notice that the terms containing the vector field ujq generate a repulsive contribution. 
To find bound states of the BN matter we calculate the energy per particle, 

(37) 



Nn + N^ Pn + Pb 

at different pjv and . It is useful to introduce instead of the relative concentration 
of antibaryons 

X=^ = fS. (38) 
Nn Pn 

Due to the charge-conjugation invariance of NN matter, its energy density must be 
symmetric under the replacement pj^ <-> p^v ■ Therefore, it is sufficient to consider only 
the values x < 1 . At fixed %, using Eqs. (27)-(35) one has in the low density hmit 

lime= "^^ + "^^^ . (39) 

PN^O 1+x 

By definition, the total binding energy of the system is equal to 

BE^{eo-e)-{NN + N^). (40) 

In the case of A^A^ matter = 'm]\f . To investigate qualitatively properties of bound 
B + ^^O systems (see Sect. IV), we study separately the case of x = 1/16 . 

Figures 1, 2 show the energies per particle of A^A" matter calculated within the TMl 
and NLZ models. The coupling constants of A" interactions are fixed by the G-parity 
transformation (see Eq. (19)). Different curves corresponds to different values of anti- 
nucleon concentration x- Dashed lines in Fig. 1 represent the results for the case x = 1/16 . 
It is seen that nucleonic matter becomes more bound after inserting a certain fraction 
of antibaryons. The maximal binding takes place at x = 1 i-e- for the "baryon-free" 
matter with net baryon density equal to zero. As one can see from Eqs. (33), (36), the 
repulsive vector contribution to the energy density vanishes in this case {ujo — 0), therefore, 
the behavior of e(pAr) is determined only by the counterbalance of scalar attraction and 

14 




0.85 ' ' 1 ' 1 ' 1 ' 1 ' ' ' 1 

1 2 3 4 5 6 

FIG. 2: Comparison of energies per particle calculated within the TMl and NLZ models for 
X = and 1/16. 

effects of Fermi motion. It is interesting to note that such baryon-free NN matter is 
fully symmetric with respect to interchanging nucleons and antinucleons. Thus, there is 
absolutely no reason that the and coupling constants would violate the G-parity 
symmetry. 

To illustrate the model dependence of the results, in Fig. 2 we compare predictions 
of the TMl and NLZ models for x = and x = 1/16. One can see that the NLZ 
parametrization predicts larger binding of the A^A^ matter at densities ~ 2po • The 
calculated parameters of bound states for x = (pure nucleonic matter), 1/16 and 1 are 
presented in Table II. By using Eq. (40) and the values of minimal energies per particle 
given in this table (for x = V^S), one can estimate the binding energy of a bound p +^^0 
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TABLE II: Characteristics of bound states of pure nucleonic matter (x = 0) as well as the NN 
matter predicted by the TMl and NLZ models. 



Model 


TMl 


NLZ 


X 





1/16 


1 





1/16 


1 


Pn/po 


0.99 


1.75 


4.99 


1.01 


1.84 


4.04 


e (MeV) 


922 


875 


449 


923 


860 


412 


(MeV) 


593 


393 


90 


544 


224 


42 



system. This leads to the result 

^ f 1070 McV (TMl), ^ ^ 

BECIO)^{ (41) 
I 1330 MeV (NLZ) . 

On the other hand, applying the formalism of Sect. II to a finite system ^^O [24] gives for 
its binding energy 1159 MeV (TMl) and 828 MeV (NLZ2). Comparison of these results 
shows influence of the finite size (surface), rearrangement (polarization) and Coulomb 
effects which are not taken into account in the infinite matter calculations. 

More detailed results for the NN matter are represented in Fig. 3 in the form of contour 
plots of e{pN,Pj[) on the plane pn — Pn calculated within the TMl model. The states 
with fixed x lie on straight lines going from the origin of the plane. The maximal binding, 
about 490 MeV per particle, is predicted for symmetric systems with pjy = p^r ~ 5po- 
Most likely, the hadronic language is not valid at such high densities and the quark- 
antiquark degrees of freedom are more appropriate in this case. 

Strong binding effects can also be expected in a deconfined qq matter. In Refs. [33, 34] 
we have performed calculations for such matter within the generalized Nambu-Jona- 
Lasinio (NJL) model. As an example. Fig. 4 shows energy per particle, e = E/{Nq + Nq), 
for nonstrange isospin-symmetric qq matter at zero temperature. Here Nq {Nq) denotes 
the total number of u, d quarks (antiquarks) . The figure displays contour plots of e {pq, pq) 
predicted by the SU(2) NJL model (for details, see Ref. [33]). The binding in this case 
is generated by an attractive interaction in the scalar-pseudoscalar channel. Again, one 
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2 4 6 8 10 

FIG. 3: Contours of energy per particle (shown in GeV near the corresponding curves) on the 
Pn — Pj; plane, calculated within the TMl model. The dashed line corresponds to states with 
X = 1/16 . Solid point shows position of minimal e on this line. 

can see that the strongest binding is predicted for the baryon-symmetric case pq = pq. 
The maximum binding energy per qq pair is 

BE ~ {rrfq^" - 270 MeV) x 2 ~ 60 MeV , (42) 

where m™^ ~ 300 MeV is the constituent mass of light quarks in the vacuum. Approx- 
imately three times larger binding energies have been found in Ref. [34] for cold quark- 
antiquark matter with admixture of strange s, s quarks. On the basis of this finding, in 
Refs. [33, 34] we have predicted possible existence of new metastable systems, mesoballs, 
consisting of many quarks and antiquarks in a common "bag". Such systems are char- 
acterized by a small (or zero) baryon number and their decay should occur mainly via 
emitting pions from the surface. It is interesting to note that multi-quark-antiquark 
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2 4 6 8 10 



FIG. 4: Contours of energy per particle, e = E/[Nq + Nq) , (shown in GeV near the correspond- 
ing curves) of cold quark-antiquark matter calculated within the SU(2) NJL model [33]. The 
horizontal (vertical) axis corresponds to the total density of u + d quarks (antiquarks) in units 
of pq = 0.15 fm~^ . 

clusters were also predicted [35] on the basis of the MIT bag model. It is instructive to 
compare the value (42) with the binding energy of baryon-symmetric NN matter cal- 
culated within the RMF model. Using Table II for the case x = 1 5 we get the binding 
energy per quark-antiquark pair, 2(mAr — e)/3 ~ 326 (350) MeV, within the TMl (NLZ) 
model. Therefore, the hadronic approach may overestimate real binding energies of cold 
and dense baryon-free matter by a factor ~ 5 . 

As noted above, the G-parity symmetry may be violated in a dense baryon-rich mat- 
ter. As a consequence, the G-parity motivated choice of antibaryon coupling constants, 
Eq. (19), may overestimate binding energies of the BN matter. To elaborate on this 
issue, we have performed an analogous calculation, but with reduced couplings g^-g , gj^ 
as defined in Eq. (20). The results of calculations within the TMl model for different 
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FIG. 5: Energy per particle of the NN matter versus the nucleon density at fixed x = 1/16. 
Different curves correspond to the TMl calculation at different values of the parameter ^ char- 
acterizing deviation of the TV coupling constants from the G-parity motivated values. 

values of ^ are presented in Fig. 5. It shows the energy per particle of NN matter with 
X = 1/16. The lower curve corresponds to the case of the exact G-parity symmetry. One 
can see that a reduction of ^ results in a smaller binding and compression of the NN 
matter. However, the admixture of antinucleons becomes relatively unimportant only at 
very small antinucleon couplings, corresponding io ^ < 0.25 . It is shown below that the 
same conclusion follows from more refined calculations for the finite ^^O system. 



TABLE III: The parameters of the TMl and NLZ models in the hyperonic sector. The remaining 
parameters are given in Table I. 





mj^ (MeV) 


m„* (MeV) 


(MeV) 


9c-K 


9ujk 


9p-K 


9a* A 


9(I>A 


NLZ 


1116 


975 


1020 


6.23 


-8.61 





6.77 


6.09 


TMl 


1116 


975 


1020 


6.21 


-8.41 





6.67 


5.95 



19 



We end this section by considering cold nucleonic matter with admixture of antihy- 
perons y = A, E . . . As proposed in Ref. [36], the observed data on AA interaction can 
be reproduced within the RMF model by including additional scalar (cr*) and vector (0) 
meson fields coupled only to hyperons. By the same reason, to take into account the 
interaction between antihyperons in the YN matter, we generalize the Lagrangian (4) by 
introducing the additional terms [69] 

+"0 Y (^a. Y - Qc^Y (t>^lp) ■ (43) 

In our calculations of the h.N matter and bound A-nuclear systems we use the values 
of parameters suggested in Ref. [36]. The A-meson couplings were obtained from the 
A-meson coupling constants by using the G-parity transformation. 



TABLE IV: Characteristics of bound states of pure nucleonic matter (% = 0) as well as the 
AA^-matter predicted by the TMl models. 



X 





1/16 


1/4 


1 


Pn/po 


0.99 


1.47 


2.49 


4.08 


eo - e (MeV) 


16 


44 


135 


252 


mlf (MeV) 


593 


456 


257 


44 


mi- (MeV) 


905 


812 


662 


420 



It is easy to see that the hyperon-hyperon interactions lead to the following effects. 
First, the gap equation for the antihyperon effective mass is modified (as compared 
to Eq. (29)) due to the additional cr* field: 

= mA - g^Y(^ - t/^^yfJ* . (44) 

Second, the energy density of the YN matter contains now the additional contribution of 
(7* and 4> mesons, 

SeYY^^imlal + mlcPl) . (45) 



20 



.Q 4 I I I I I I I I I I 

■ 2 4 6 8 10 

FIG. 6: Energy per particle of AA^ matter versus the nucleon density calculated within the TMl 
model. Different curves correspond to different values of the A concentration x = Pa"/ Pn ■ The 

curves are shifted along the vertical axis by a ^"dependent value eo = lim e . Point A marks 

PJV^O 

the position of energy minimum for x = 1/16 . 

Figure 6 shows the energy per particle of the AA^ matter, calculated within the TMl 
model, generalized in accordance with Eqs. (44)-(45), with parameters listed in Ta- 
bles I, III. To compare results for different x ? the energy per particle is shifted by a 
corresponding vacuum value eo defined in Eq. (39). The calculated parameters of bound 
states are given in Table IV. Comparison with results obtained earlier for the NN matter 
shows that binding energies and equilibrium densities are noticeably smaller in the AA^ 
case. For example, the binding energy predicted for the ^.^O "nucleus" equals approx- 
imately 44 X 17 ~ 748 MeV which is about 30% smaller than for ^^O bound state. 
This difference is explained by a smaller scalar coupling of A particles as compared to 
antinucleons {g^-x/g^N - 2/3) . 
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IV. FINITE ANTIBARYON NUCLEAR SYSTEMS 



A. Light nuclear systems with antiprotons 

The nucleus ^^O is the lightest nuclear system for which the RMF approach is con- 
sidered to be reliable. This nucleus is included into the fit of the effective forces NL3 
and NLZ2. Therefore, we choose this nucleus as a basic system to study the antibaryon- 
nuclear bound states. First, let us consider an antiproton bound in a ^^O nucleus. For 
clarity we use the notation ^^O for such a combined system. A priori it is unclear which 
quantum numbers has the lowest bound state. We first assume that this is the state 
and later on we shall check this assumption. Results of our self-consistent calculations 
for both ^^O and are presented in Fig. 7. It shows 3D plots of the nucleon densities 




FIG. 7: 3D plots of nucleon density in the ^^O nucleus (left) and in the bound system 
(right) calculated within the NL3 model. 

for the G-parity motivated antiproton couplings (^=1). One can see that inserting an 
antiproton into the nucleus gives rise to a dramatic rearrangement of nuclear structure. 
This effect has a simple origin. As explained above, the antiproton contributes with the 
same sign as nucleons to the scalar density (see Eq. (21)), but with the negative sign to the 
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vector density. This leads to an overall increase of attraction and decrease of repulsion for 
surrounding nucleons. To maximize attraction, protons and neutrons move to the center 
of the nucleus, where the antiproton has its largest occupation probability. This leads to 
a strong compression of the nucleus. 



0.6 




NLZ2 

--- -- NL3 


So.4 




TMl 


0.2 






0.0 




' 1 ' 1 ' 1 ' 1 ' 1 ' 1^ — 






0.6 




/ 


0.4 




' - \ 
/.-'■-, 
' '■ 1 
' ■ '.1 


2 




/ \ 


0.2 




/ \ 


0.0 






-4 


-2 2 4 



r(fm) 




-4 -2 2 4 
r (fm) r (fm) 



FIG. 8: The left panel represents the sum of proton and neutron densities as function of nuclear 
radius for ^^O without (top) and with an antiproton. The left and right parts of the upper 
middle panel show separately the proton and neutron densities, the lower part of this panel 
displays the antiproton density (with minus sign). The right panel shows the scalar (negative) 
and vector (positive) parts of the nucleon potential. Small contributions shown in the lower row 
correspond to the isovector (p-meson) part. 

Figure 8 shows the densities and potentials for ^^O with and without the antiproton. 
For normal ^^O all RMF parametrizations considered produce very similar results. The 
presence of the antiproton changes drastically the structure of the nucleus. Depending on 
the parametrization, the sum of proton and neutron densities reaches a maximum value 
of (2 — 4) po ) where po ~ 0.15 fm~^ is the normal nuclear density. The largest compression 
is predicted by the TMl model. This follows from the fact that this parametrization gives 
the softest equation of state as compared to other forces considered here. According to 
our calculations, the difference between proton and neutron densities is quite large, which 
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leads to an increase in symmetry energy. The reason is that protons, though they feel 
additional Coulomb attraction to the antiproton, repel each other. As a consequence, 
neutrons are concentrated closer to the center than protons and the symmetry energy 
increases. 
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FIG. 9: Left panels: radial profiles of the vector (pj) and scalar {psj) densities of antipro- 
ton (j = p) and nucleons (j = N) in the ground state of ^^O as predicted by the NLZ2 and NL3 
calculations. Right panels show the corresponding profiles of the effective nucleon mass. 

Figure 9 shows radial profiles of scalar and vector densities of nucleons and the an- 
tiproton in ^^O . As compared to the normal ^^O nucleus, absolute values of vector and 
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scalar densities increase in the central region of the nucleus. This leads to a strong drop 
of the effective nucleon mass near the nuclear center [70], which in turn suppresses the 
local annihilation rate of the antiproton (see Sect. V). 




FIG. 10: Proton (left), neutron (middle) and antiproton (right) energy levels for the nucleus 
^^O with one antiproton and without it (rightmost columns in the left and middle panels). 

Since nucleons feel a deeper potential due to the presence of the antiproton, their 
binding energy increases too. This can be seen in Fig. 10. The nucleon binding is largest 
within the NL3 parametrization. In the TMl case, the sl/2+ state is also deep, but 
higher levels are less bound as compared to the NL3 and NLZ2 calculations. This is a 
consequence of the smaller spatial extension of the potential in this case. The highest 
sl/2~ level in the TMl calculation (see Fig. 10) is even less bound than for the system 
without an antiproton. 

For the antiproton levels, the TMl parametrization predicts the deepest bound state 
with binding energy of about 1130 MeV. The NL3 calculation gives nearly the same 
binding, while in the NLZ2 case, antiproton levels are more shallow and have smaller 
spacing. It should be noted that the antiprotons are more strongly bound than was 
obtained in Ref. [15]. This follows from the fact that here we consider both nucleons 
and the antinucleon self-consistently allowing the target nucleus to change its shape and 
structure due to the presence of the antiproton. The total binding energy of the ^^O 



25 



system is predicted to be 828 MeV for NLZ2, 1051 MeV for NL3, and 1159 MeV for 
TMl. For comparison, the binding energy of the normal ^^O nucleus is 127.8, 128.7 and 
130.3 MeV in NLZ2, NL3, and TMl, respectively. 




FIG. 11: 3D plots of nucleoli densities for the nucleus ^Be without (left) and with (right) 
antiproton calculated within the NL3 model. 

As a second example, we investigate the effect of a single antiproton inserted into 
the ^Be nucleus. In this case only the NL3 parametrization was used (the effect is similar 
for all three RMF models). The normal *Be nucleus is not spherical, exhibiting a clearly 
visible a — a structure with the deformation — 1.20 in the ground state. As one can see 
from Fig. 11, inserting an antiproton gives rise to the compression and change of nuclear 
shape, which results in a much less elongated nucleus with ^2 — 0.23. Its maximum 
density increases by a factor of three from 1.3po to 4.1po • The cluster structure of the 
ground state completely vanishes. A similar effect has been predicted in Ref. [37] for 
the K~ bound state in the ^Be nucleus. In our case the binding energy increases from 
52.9 MeV (the experimental value is 56.5 MeV) to about 700 MeV. 
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FIG. 12: The profiles of net baryon density, Pn — Pp, in the p+^He (left) and (right) 
systems, calculated within the NLZ2 model. 

Figure 12 represents profiles of the net baryon density i.e. the sum of the proton and 
neutron densities minus the antiproton density. For both considered systems, there is a 
dip in the center surrounded by the region with baryon density increased as compared to 
normal nuclei. 



TABLE V: Total binding energies of p+ O system with antiproton occupying different states. 



TTT 


1/2+ 


3/2- 


1/2- 


1/2- 


NL3 


-1051 


-1008 


-920 


-780 


NLZ2 


-828 


-938 


-804 


-876 



Now we address the question of wether the antiproton 1/2+ state is energetically the 
most favorable configuration of the system or not. To answer this questions, we put 
the antiproton into various states and calculate the total binding energy. The results 
are shown in Table V. In the case of NL3, the lowest 1/2+ state indeed corresponds to 
the configuration with lowest energy. However, within the NLZ2 calculation a state with 
the antiproton in the 3/2~ state corresponds to the ground state of the system. The 
3/2" state has a different spatial distribution and thus, leads to a different shape of the 
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FIG. 13: 3D plot of nucleon density in the in the bound A + system calculated with the NLZ 
model. 

potential felt by nucleons. So, even though the antiproton in the 3/2^ state is less bound 
than in the 1/2"'" state, nucleons become more bound in this case and the binding energy 
of the whole system increases. 

B. Light systems with antilambdas 

Bound states of nuclei with antihyperons are especially interesting because we expect 
longer life times in this case (see Sect. V). We performed calculations of the He and ^.^O 
systems within the NLZ model. The assumed A couplings with the meson fields are given 
in Table III. 

Figure 13 shows the 3D plot of the sum of proton and neutron densities in the ^^O 
nucleus. Comparison with the results for the ^^O systems (see Fig. 7) shows noticeably 
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FIG. 14: Radial profiles of the net baryon density, — Pj^-, in the A -|-^He (left) and A -|-^^0 
(right) systems calculated with the NLZ model. 

lower nucleon densities in this case. However, we still predict significant compression of 
nuclei containing antihyperons. Figure 14 gives profiles of the net baryon density in the 
i-Hc and ^-^O nuclei. It is interesting that the central density dip is even more pronounced 
here as compared to nuclei with antiprotons (see Fig. 12). This effect can be explained 
by a stronger localization of the antilambda due to its larger effective mass. 



C. Calculations with reduced antibaryon couplings 

As was pointed out earlier, the G-parity transformation may not work on the mean- 
field level. Therefore, we have performed calculations of bound p-\-A systems with reduced 
antiproton couplings to mean meson fields. Figure 15 shows the results for the ^^O nucleus 
obtained for several values of the parameter ^ introduced in Eq. (20). 

One can see that at ^ ^ 0.25, the compressed shape of the nucleus is not very much 
affected. The largest compression, of course, occurs for ^ = 1, i.e. in the case of the 
exact G-parity. At ^ smaller than 0.05 (NL3) or 0.01 (NLZ2), the nucleon density profiles 
practically coincide with the density distribution of the normal ^^O nucleus. On the 
other hand, the antiproton density profile becomes rather flat in this case. Somewhere in 
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FIG. 15: The profiles of densities and mean-field potentials in the bound system calcu- 

lated with reduced antiproton coupling constants. The upper (lower) panel corresponds to the 
NL3 (NLZ2) model. Different curves correspond to different values of the parameter ^ defined 
in Eq. (20) . Shown are the nucleon densities (left) , antiproton densities (middle) and scalar and 
vector potentials for nucleons (right). Note, that at ^ < 1 the antiproton potentials are smaller 
due to reduced p couplings. 

between, there exists a critical ^ value separating these two regimes. Detailed calculations 
with different ^ show that there is no smooth transition but rather an abrupt jump. This 
is an indication that, to a large extent, the shell effects control the structure of these 

systems. 

The results of these calculations show that the conclusions made in Ref. [24] and 
in the present work do not strongly depend on the actual values of g^jf and g^jf. As 
a consistency check, we verified that at ^ — > the binding energy of ^^O goes to the value 
of normal ^^O [71]. This is demonstrated in Fig. 16 where one can see strong increase of 
the binding energy for ^ > 0.1 . 

Similar results take place for the p+^He system. They are presented in Fig. 17. Of 
course, this is system is rather small and the mean-field approximation may be not ac- 
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FIG. 16: Binding energy of ^^O calculated in the TMl and NLZ models as a function of the 
parameter ^ . 
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FIG. 17: Same as in Fig. 15 but for ^He with an antiproton. 
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curate in this case. Nevertheless, we expect to get a quahtative picture even for such a 
hght system. At ^ = 1 the nucleon density in the pie nucleus reaches a maximum value 
of about 6po which is noticeably larger than in ^^O. As in the case of ^^O, depending on 
the ^ value, there are two different configurations of the bound system, namely a highly 
compressed one and the other resembling normal helium plus a quasi free antiproton. For 
both RMF models, the reduction of antiproton couphng constants up to ^ ~ 0.25 does 
not produce a strong effect in the density distribution. 



D. Rearrangement of nuclear structure and Dirac sea 




FIG. 18: Proton and antiproton potentials as well as corresponding single-particle levels for 
the systems ^He, p+^Re (left side) and ^^O, (right side) calculated within the NLZ2 

model. The full and open dots in the Dirac sea show the most bound occupied and vacant states, 
respectively. 

The presence of an antiproton within a light nucleus leads to a drastic rearrangement 
of its structure and to the polarization of the Dirac sea. We illustrate this effect in 
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Fig. 18, where the effective potentials and single-particle levels are shown for helium 
and oxygen, in each case with and without an antiproton. A hole at the deepest bound 
level in the lower well corresponds to a positive energy antiproton. After inserting an 
antiproton the nuclear structure is strongly rearranged. Both the nucleonic as well as the 
antiparticle states become deeper which gives rise to increase of the binding energy by 
several hundred MeV. The lower well in the case of the p+^He system is very narrow 
and its bottom overlaps with the positive energy well. The uncertainty relation prevents, 
however, the strongest bound antiparticle and particle states to be very close in energy. 
The lower well in the case of oxygen also becomes narrower as compared with normal ^^O 
nucleus. 

The rearrangement of the Dirac sea due to the presence of an antiproton leads to 
increasing number of negative energy states. The spin-orbit potential for antiparticles is 
very small, because, in contrast to nucleons, scalar and vector potentials for antiprotons 
nearly cancel each other [72] . Still, there are gaps between the antiparticle levels of the 
order of a hundred MeV (see Fig. 10). This is a consequence of a smaller width of the 
lower well as compared to the negative energy well in the normal nucleus. 

E. Heavy nuclei containing antibaryons 

In this section we investigate the effect of antiprotons and antilambdas inserted into 
heavy nuclei. Because of larger radii of these objects {R^7 fm), one may expect appear- 
ance of a local compression zone in the central region instead of a more homogeneously 
compression in lighter systems. This is exactly what follows from our calculation. Fi- 
gure 19 shows the sum of proton and neutron densities, for the case of doubly-magic 
lead nucleus with one deeply-bound antiproton. Due to the presence of an antiproton, 
a small core of highly compressed nuclear matter appears at in the center of the ^^^Pb 
nucleus. In addition, the lead nucleus becomes deformed and acquires a prolate shape. 
As shown in Fig. 20 a similar structural change occurs when implementing the A particle 
into the ^''^Pb target. Figure 21 shows single-particle spectra of protons, neutrons and 
antiprotons in the normal lead as well as in the bound ^°®Pb and pb nuclei. One can 
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FIG. 19: 3D plots of nucleoli density in the ^ospj-, nucleus (left) and in the bound p+'^^^Ph 
system calculated within the NL3 model. 

see that implementing antiparticles results in a strong change of shell structure. Due to 
the axial deformation, the system looses its degeneracy and the shell structure becomes 
partly washed out. Only the 1/2+ and 3/2— levels exhibit large binding. The deepest 
antibaryon states have binding energies of more than 900 (600) MeV in the case of the 
^°^Pb (^5?Pb) system. On the other hand, some of the single-particle states become even 
less bound in the lead nucleus with antibaryons. 

These results can be qualitatively understood from the fact that the nucleon potential 
in the center of the bound system is very deep but narrow. Two effects prevent deeper 
binding of single-particle levels. First, a narrow potential gives rise to large uncertainties 
in particle momenta. This in turn increases the kinetic energies of particles and therefore, 
reduces their binding. Second, the single-particle states with larger angular momenta are 
mainly localized at larger radii and thus do not have much overlap with the deep central 
region of the potential. The total binding energy of the ^^^Pb nucleus predicted by the 
NL3 calculation equals 2412 MeV or 11.5 MeV per particle. This is significantly smaller 
than 61.8 MeV per particle in the ^|0 system. Furthermore, the total energy gain in the 
binding energy after inserting an antibaryon into the ^°^Pb nucleus (780 MeV) is smaller 
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FIG. 20: 3D plots of nucleoli density in the bound A+ Pb system calculated within the NLZ 
model. 
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FIG. 21: Single-particle proton (left), neutron (middle) and antiproton spectra for the systems 
2°^Pb, ^+208pb (NL3 forces) and A+^o^Pb (NLZ2 set). 
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than in the ^^O case (920 MeV). 
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FIG. 22: The profiles of nucleoli density in the Pb system at different values of the para- 
mater ^ . The left and right panels correspond, respectively, to the equatorial and coaxial planes. 

Figure 22 shows the nucleon density profiles in the Pb system calculated with 
^ = 1,0.5 and (the normal Pb nucleus). Because of the axial deformation we show 
separately the profiles for radial (equatorial) and coaxial planes. This figure demonstrates 
a very interesting behavior. Contrary to a naive expectation, reducing from 1 to 0.5 
leads to increasing density of the central core [73]. We believe that this can be explained 
by the competition of two opposite trends taking place at decreasing ^ . The first one is the 
reduction of the binding potential while the second one is the increase of the antiproton 
effective mass. The second effect leads to a stronger locahzation of the antiproton wave 
function near the minimum of the effective potential. As a consequence, the potential 
acting on surrounding nucleons increases in the central regions. 
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F. Systems containing several antibaryons 

It is amazing to consider nuclear systems with more than one trapped antibaryon. For 
instance let us consider doubly-magic oxygen containing different numbers of antiprotons, 
namely Np = 2,4,6,8 and 10. The corresponding baryon densities are shown in Fig. 23. 



0.3 
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FIG. 23: Net baryon density, pN — Pp, for the ^^O nucleus with 2, 4, 6, 8 and 10 antiprotons, 
calculated within the NL3 model. 

The total binding energies for the systems containing 2, 4, 6, 8 and 10 antiprotons are 
1541, 2792, 3847, 5006, and 6300 MeV, respectively. As we can expect, with increasing 
number of antiprotons, the region of reduced net baryon density becomes broader. In the 
case Np = 10 one can see a qualitative change of nuclear structure, with appearance of 
a baryon density peak in the center region and a zone of reduced density at the nuclear 
periphery. 

One can also think of most extreme case of finite systems with equal numbers of baryons 
and antibaryons, i.e. about systems with total baryon number B — . These systems can 
be called self-conjugate nuclei, since their charge conjugation leads to the same object. 
Let us first consider the bound a — a system. As seen in Fig. 24, the (anti) particle 
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FIG. 24: The system a— a calculated within the NLZ2 (solid line) and NL3 (dotted line) models. 
Shown are baryon and antibaryon densities (left) as well as the scalar potential and energies 
of the lowest singic-particlc states (right). The dashed line shows the nucleon density of the 
normal a particle calculated within the NLZ2 model. 

densities in this system reach about lOpo • Of course, one should expect a breakdown of 
the RMF model with nucleonic degrees of freedom at such high densities. Most probably 
the (anti)baryons will dissolve into a deconfined state of cold and dense qq matter (see 
Sect. III). 

Due to the full symmetry between baryons and antibaryons (or quark and antiquarks) 
such systems have both vanishing baryon and electric charges. Accordingly, the vector 
fields, namely the lo meson and Coulomb fields, vanish too. Only the scalar density and 
with it the scalar potential are nonzero in a symmetrical system. Since nucleons and antin- 
ucleons have the same scalar potential, they can occupy the same single-particle states. 
In this case the spatial parts of the wave functions of protons, neutrons, antiprotons and 
antineutrons are identical. Hence, the spatial overlap between particles and antiparticles 
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FIG. 25: The system ^^0-^*^0 calculated within the NL3 model. The dotted lines show the 
baryon and antibaryon densities (left) as well as the scalar potential and the lowest single-particle 
states (right). The dashed line represents the nucleon density of the normal ^^O nucleus. 

is maximal. On the other hand, due to the very high density of baryons and antibaryons, 
the effective masses are very small in the core region. The calculated binding energy of 
the a - a system is very large, it is 2649 MeV for NL3 and 2235 MeV for NLZ2. This 
is about 100 times the value for a single a. particle and 40 times the binding energy of 
normal ^Be nucleus. 

As a second example we consider the ^^O-^^O system. The calculated profiles of 
(anti)nucleon density and scalar potential are shown in Fig. 25 [74]. In this case we predict 
formation of a highly compressed system with maximal nucleon and antinucleon densities 
of about 5po • One can see in Figs. 24-25 that radial profiles of the scalar potential in 
symmetric A^A systems are noticeably broader as compared to the (anti)nucleon density 
distributions. Such a behavior follows from the finite range of a meson field (see Eq. (11)). 
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V. ANNIHILATION OF DEEPLY BOUND ANTIBARYONS 



A. NN annihilation in vacuum 

At low energies the total cross section of pp annihilation in vacuum can be well ap- 
proximated as [39] 

= ^ + V- ' 

where v^ei is the relative velocity of p with respect to the proton, C = 38 mb, D — 35 mb-c . 
In the following we neglect the isotopic and Coulomb effects assuming that pp and pn 
annihilation cross sections are approximately equal. Using this approximation we discuss 
below the isotopically averaged A^A^ interactions. 

When the NN cm. energy ^/s is close to the vacuum threshold value, i.e. at 
2m N — 1.88 GeV, the relative velocity becomes small and the NN annihilation 
cross section is approximately proportional to 1/i'rei [75]. This limiting case of annihila- 
tion "at rest" {v^-ei 0) has been extensively studied in Refs. [41-47]. Experimental data 
on many exclusive annihilation channels pp ^ c ^ n^7r are also available. Here interme- 
diate states c include direct pions as well as heavy mesons r7(547), p(770), c<j (782) . . . 

Various channels of the pp annihilation at rest are listed in Table VI. They are sorted 
into groups with different final pion multiplicities . The most important exclusive 
channels c are given in the second column. In the third column we give the corresponding 
threshold energies in the cm. frame (i.e. the total mass of particles in the intermediate 
state). The fourth column shows the observed branching ratio of a given channel c: 

= <7,,_,/a|7 . (47) 

In the last column we give references to publications where experimental information on 
branching ratios has been found. The subscripts "dir" in the second column mean that 
only nonresonant contributions are given in the corresponding line. For example, the 
branching ratio of the (7r"^7r~7r'')dir channel has been calculated using the formula 

-Spp— »(7r+7r~7rO) -^pp— ^Tr+TT-Tr" -^pp— »7r0p0 -^pp— »7r±pT j (^8) 
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TABLE VI: Exclusive channels of pp annihilation at rest in vacuum 





c 


V^,hr(GeV) 


Sc(%) 


Refs. 


z 


9-7r0 


u.z < 


u.u t 






/i /i 


f) 28 


n 31 




Q 
O 




n /II 


n 7 


fzL7l 




^TT TT TT jdir 




l.o 


[40J 




vr p 


n Qi 


1 7 


[40J 






n Qi 
u.y i 




[40J 


4 


(4^°)dir 


0.54 


0.5* 


[47] 




(7r+7r-27r0)di,. 


0.55 


7.8'= 


[41, 42, 46] 




^/vr /vr jdir 


u.oo 


^.Z 


[4Z, 4dJ 




71 W 


n Q9 


U.D 


[4bJ 




/I /I ^ 


l.U(J 


O.L) 


[4ZJ 








n Q 
u.y 


[4ZJ 


5 


(57rO)dir 


0.68 


0.5^ 


[44, 45] 




(7r+7r"37r°)dir 


0.69 


20.1-'? 


[41, 44, 47] 




(27r+27r-7r0)dir 


0.70 


10.4'' 


[42, 44, 46] 




27r°r? 


0.82 


0.7 


[44] 




TT TT // 


u.oo 


i.O 


[44J 




ZTT u/ 


l.UO 


9 fi 
Z.O 


f47l 
[4CJ 






l.UD 


0.0 


[40J 




fj uJ 




9 


[40J 


6 


(7r+7r-47r0)dir 


0.82 


1.9' 


[41, 42] 




^ZTT ZTT ZTT J dir 


u.oo 


io.o 


\A^ 471 




\6'K OTT j dir 


u.o^t 


9 n 
z.u 


f49l 
[4ZJ 






1.32 


1.5 


[47] 




2a; 


1.54 


3.0 


[47] 


7 


(27r+27r-37r0)dir 


0.97 


4.0^= 


[41, 47] 




(37r+37r-7rO)dir 


0.98 


1.9 


[42] 






1.47 


1.0 


[47] 



°- see Eq. (48) 

average between two values given in Ref. [47] 

calculated as Bpp->7r+7r-27r'' ~ Bpp^^^o^ — Bpp^p+p- 

calculated as Bpp^2n+i-K- — Bpp^.^+.^-pO 
^ average between two values referred in Ref. [42] 

calculated as Bpp^^^a — Bpp_^2TT'>riB n—^sn^ 
^ calculated as Bpp^^+^-^^o — i3pp^27i-Or)-B^^^+^-^o — Bpp^^+^-^Bn^^T^o — Bpp^^TT^u 

^ calculated as Bpp_y2Tr+2TT-Tr° ~ -^pp— vTr+Tr-T?-^??— ►Tr+Tr-Tr" ~ -Spp^TT+Tr-oi ~ -^pp— >p°a; 

' calculated as Bpp^^+^-^^o („>2) - Bpp^„+„-2^o - Bpp^^+^-^^o using data from Ref. [42] 

^ calculated as -Bpp_>27r+27r-27r° ~ -Bpp^wjy-Sjj^TT+Tr-TrO ~ -Bpp^2w 



where the first term in the r.h.s. is the total branching ratio of the reaction pp — > 7r"^7r~7r° 
and the two other terms give the contributions of the up intermediate states. Here we 
take into account that p — > 27r is the dominant channel of the p meson decay [48] . The 
channels c = 37r°, 47r°, tt'^tt'^ti^, Stt+Stt" and 37r"'"37r~7r° are classified as direct, because 
no resonance contributions have been found there. In Table VI we include only main 
annihilation channels with Be > 0.5% (for > 2). We do not separate resonances in 
channels with total number of intermediate mesons exceeding three. 

B. Kinetic approach to annihilation in medium 

Let us first consider annihilation of slow antibaryons, B , in homogeneous nucleonic 
matter with a small admixture of antibaryons. Using again the mean-field approach we 
treat the B + N matter as a mixture of ideal Fermi gases of nucleons and antibaryons 
interacting with mean meson fields. Assuming that annihilation on multinucleon clusters 
is relatively small (see below), we consider here only the binary BN annihilation. Then 
the local rate of B annihilation per unit vohimc can be written as a sum over all BN 
annihilation channels with production of n > 2 mesons in the final state: 

dtdV ~ ^ dtdV ■ ^ ' 

n 

Here Mi{l = 1, . . . ,n) denote secondary mesons. In the following we disregard in-medium 
modifications of these mesons, assuming that their 4-momenta ki satisfy the mass-shell 
constraints, kf — mf , with vacuum masses mi [76]. 

In the quasiclassical approximation one can describe the phase-space density of an- 
tibaryons {i — B) and nucleons {i — N) by distribution functions /j(x, p) , where x — {t, r) 
is the space-time coordinate and p is the 3-momentum of the ^-th particle. Within the 
mean-field model, the kinematic part of the single-particle energy is written as 

Ei = ^mf + p', (50) 

where m* is the effective mass of the i-th particles {i — B,N). Their vector and scalar 



42 



densities are determined by integrals over 3-momenta: 

Pi{x)^ I d'pMx,p), (51) 



/Tn* 
d'p^J^{x,p). (52) 

At zero temperature we use Fermi distributions fi — Q(pFi — \p\) , where is the 

spin-isospin degeneracy factor of the i-th particles, ppi — {^'k'^ pi/viY^^ is their Fermi 
momentum and Q{x) = (1 + sgna;)/2 [77]. 

Within the kinetic approach (see e.g. Ref. [52]), the rate of the reaction 
BN c {c = Ml . . . Mn) can be written as 



dtdV 
^d% 



f d-^PB f ( ^ f d-^PN f ( V 

J -^fBi^^PB) J -^fN{x,PN) X 



\l^W,{P^p^ + Pr,\ki,...,K)5^^\P-ki-...-K). (53) 

1=1 ' 

Here Wc is the BN — > c transition probability, pi — {Ei,Pi) and ki — {Ei,ki) are the 
4-momenta of the initial particles and the final mesons. 

Following the standard arguments of a statistical approach we assume that the transi- 
tion probabilities do not depend sensitively on the particles' momenta so that W^. can be 
replaced by constants, 

W, (P| . . . , k^) c±W,^ const . (54) 

Then the contributions of different reaction channels are determined mainly by their 
invariant phase space volumes 

\W-^^''\P-Y.^^)^ (55) 

J 1 1 ^ 7 1 



1 = 1 ' 1 = 1 



where s = P = {p^ + pjv) is the cm. energy squared available for the reaction. 
Now Eq. (53) can be written as 

where Fc is the partial annihilation width for channel c. The angular brackets denote 
averaging over the momentum distribution of incoming particles, 
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where O is an arbitrary function of 3-momenta Pb,Pn ■ 

Let us consider first the hmit of dilute matter, i.e. p^i Pn 0. In this case, ppi —>■ 0, 
m* rrii {i = B,N) and the distribution functions fi{x,p) can be formally replaced 
by Pi 5^^^ (p) . As a consequence, at p^, ^ we get the following relation: 

5?^.^,,^-.^. (58) 
ataV rrv^rriN 

Here Sq = {m-g + mNY corresponds to the energy available for BN annihilation at rest in 
vacuum. On the other hand, in this limit the reaction rate can be expressed [52] through 
the vacuum cross section of the BN annihilation at Vrei 0: 

^^^y - PbPN [CTBN^c'^rel) , (59) 

where the subscript " 0" indicates that the quantity in brackets is taken at Pb^n ~^ ■ 
By comparing Eqs. (58) and (59) one can express the transition probabilities Wc by 
experimental cross sections of the BN annihilation: 

In a general case of finite densities, the partial annihilation width, defined in Eq. (56), 
can be expressed as 

rc = r„.(A.w!^), (61) 

where 

Toc = PN (c^SAT-^c^rel) , (62) 

is the axillary width, corresponding to the vacuum cross section, and 

- li («^' 

is the phase-space suppression factor, which plays a central role in our estimates. 

In nuclear medium the single-particle energies of antibaryons and nucleons are modified 
by the scalar and vector potentials. The minimum energy oi a, BN pair available for 
annihilation, i.e. the reaction Q value, is 

Q = ^/s; = + V-B + m% + Vn , (64) 
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where m* and Vi are effective masses and vector potentials of particles {i — B,N). In the 
case of antinucleons with the G-parity transformed potentials = m*j^ and Vjf = — V/v 
hence Q = 2m^ . 

This value can be significantly reduced as compared to the minimal energy in vacuum, 
■y/so , which leads to strong suppression of the available phase space for annihilation prod- 
ucts. Actual cm. energies of BN pairs have a certain spread due to the Fermi motion of 
nucleons. It easy to show that even at ~ 2po the variation of the Q values does not 
exceed 10%. Taking into account that uncertainties in scalar and vector potentials are of 
the same order, we simplify the general expression (61) by replacing Ac(s) with Ac(s*) . 
Then the partial annihilation width in the medium can be represented in a simple form 

Tc ~ rocAc(s*)'%7v (65) 

Here 

J^^^<h{x)M^), U^) = (^) = ^P^ {i^B,N). (66) 

\Ei/ Pi 

Returning to the initial expression (56) we see that the reaction rate is proportional to 
the product of the scalar densities of annihilating particles. In fact, this is required by 
the Lorentz invariance of Eq. (56), because factors Wc and Ac are Lorentz invariants. 

The direct calculation of the scalar densities psi using Fermi distributions in Eq. (52) 
yields 

= (67) 

2 PFi PFi 

where ^{z) is the dimensionless function defined in Eq. (35). For systems with small 
admixture of antibaryons, i.e. when pp-^ <^ pFi\f,m^, one can use the approximate 
formulae 

Pfb 3 

Here is determined for pure nucleonic matter with p^ — . 



C. Phase space suppression factors 



Let us now consider in more detail the phase space suppression factors Ac introduced 
by Eq. (63). In the particular case of the two-body annihilation channels Ac can be 
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FIG. 26: Phase space suppression factors for NN annihilation into two mesons as functions of 
cm. energy in the subthreshold region. The vertical dashed lines indicate ^/s^ values adopted 
in the calculations for infinite nuclear matter (see subsection D). 



calculated analytically. Direct calculation using Eq. (55) gives 

D f \ o \^ (^1 + ^2)^ ] L (mi-ms)^] ^ 

%iv^MiM2(^) = 27rW 1 1 e(Vs-mi-m2), (69) 



where mi and m2 are meson masses. Figure 26 presents the results for the two-body 
channels involving pions and the heavy mesons r), p,u! . One can see a strong decrease of 
Ac(s) in the subthreshold region, s < sq, which is getting more prominent for heavier final 

states. 

For channels with more than two secondary mesons, the multidimensional integrals in 
Eq. (55) are evaluated numerically using the Monte Carlo method. Figure 27 shows the 
energy dependence of factors Ac for nonresonant channels NN {n'K)diT with different 
pion multiplicities n . One can see that at given s the reduction of phase space becomes 
more and more important with increasing n . In the case of the lightest annihilation 
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FIG. 27: Same as Fig. 26, but for NN annihilation into n direct pions. Note the change to a 
log scale on the vertical axis. 

channel, c = 27r, the phase space is only slightly reduced in the subthreshold energy 
region. The noticeable deviation of X2tt from unity takes place only in the vicinity of 
the 2n threshold, i.e. at ^/s —>■ 2m.j^ . On the other hand, for n > 2 factors Xmz^s) 
decrease very rapidly with decreasing s . For example, at ^/s — 1 GeV, the phase space 
factors A„7r — 0.97, 0.19, 0.014 for n = 2, 3, 4 , respectively. 

It is interesting to note that the trend in behavior of the phase space factors above the 
threshold, ^/s > 2miv, changes to opposite, i.e. the multi-pion channels become more and 
more important with increasing ^/s . This observation was used earlier [53] to explain fast 
chemical equilibration in nuclear collisions. Thus, these two trends are complementary, 
but not inconsistent with each other. 
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D. Annihilation widths in nuclear medium 



To calculate the factors Fqc which enter into the in-medium annihilation widths (65) 
we use vacuum branching ratios, = ctbn^c/ '^^n " antinucleons {B — N) 

one has 

Toe = B,pj,{a^^v,,{)^ = ToB,^ MeV . (70) 

Po 

The numerical value Fq = 104 MeV is obtained when the parametrization (46) is used for 
the total A^A^ annihilation cross section. In the calculations below we use the experimental 
Be values from Table VI. 

As an illustration, let us consider the NN — > c partial widths for the following two 
annihilation channels: c — 2t: and c — irp . The in-medium widths are calculated from 




FIG. 28: The NN c {c = 27r, irp) partial widths in nucleonic matter with density pN ■ 

Eqs. (65)-(70), in the limit pj^ ^ . The corresponding branching rations are obtained 
from the obvious relations 527r — B2t^o + B^+^- and B^^p — B^OpO + B^±p^: , using data from 
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Table VI. To study the model dependence, the predictions of two RMF models, NLZ [27] 
and TMl [29], have been compared. In both cases, the antibaryon couplings are chosen 
according to G-parity transformation, which leads to equal effective masses of nucleons 
and antinucleons. 

The widths and F^p as functions of nuclear density are shown in Fig. 28. At low nu- 
cleon densities, these widths deviate only slightly from the linear dependence Fqc oc . 
However, at Pn ^ Po, when the effective mass drops significantly, the density de- 
pendence of the in-medium annihilation widths becomes strongly nonlinear. Since the 

phase-space factors Ac(s*) vanish at ^/s^ = 2m^ < mi + m2 , both channels become for- 
bidden at large enough densities. Figure 29 shows the density dependence of the nucleon 




FIG. 29: Effective nucleon mass as a function of density predicted by the NLZ and TMl 
models. Thin horizontal lines show threshold values of for NN annihilation into 27r and np 
mesons. 

effective mass ■m%j{pN) for the same RMF models as in Fig. 28. One can see that the 
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strong model dependence of r27r is explained by a slower decrease of within the TMl 
model. It is interesting to note that all strong annihilation channels would be closed 
when m*j^{p^) < m^. In this case only electromagnetic (c = 2j) or multi-nucleon (see 
subsection F) annihilation channels would limit the life time of an antinucleon in nuclear 
medium. 

It is instructive to calculate the annihilation widths for the following three physically 
interesting cases: 

SO: Pn^Po, mlf^niN, (71) 
SI: Pn = Po, m^ = 0.65GeV, (72) 
SII : pN = 2po , m*N = 0.30 GeV . (73) 

In the first parameter set, all in-medium effects are disregarded, i.e. s* = sq, Ac = 1, 
Jjfj^ ~ 1 and, therefore, Fc ~ Fqc • In the SI case we choose the density and effective 
mass which are commonly accepted for equilibrium nuclear matter (without antinucleons) . 
Similar values are predicted by most RMF models. This case is interesting for estimating 
the annihilation width in a situation when the rearrangement effects due to the presence 
of antibaryons are disregarded [78]. 

For the case SII we take typical values for p^ and as predicted by our calculations 
for the bound ^^O nucleus. Of course, the main contribution to annihilation comes from 
the central part of the nucleus where the antiproton is localized. As can be seen in Fig. 9, 
the antiproton wave function in the lowest bound state of ^^O practically vanishes at 
r ~ it!p = 1.5 fm. The values of density and effective mass in the SII case correspond 
approximately to the values obtained by averaging the p^ and radial profiles in the 
interval r < Rp. 

Table VII presents numerical values of the A^A^ c widths Fc calculated for the above 
parameter sets using Eqs. (65), (70). In these calculations wc use the values Jjj^ ~ 2.0 (SI) 
and 7.6 (SII), which follow from Eqs. (66)-(68). In Table VII, annihilation channels are 
grouped according to the multiplicity of pions in the final state n . As compared to the 
vacuum extrapolated values (SO), the in-medium annihilation widths are significantly 
reduced for the channels with n exceeding 3 ( SI) and 2 (SII). Many annihilation channels 
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TABLE VII: Charactoristics of N annihilation in cold nuclear matter. 



Parameter set 


so 


SI 


SII 


c 


Sc(%) 


To (MeV) 


Ac 


Tc (MeV) 


Ac 


Tc (MeV) 


27r 


0.38 


0.39 


0.99 


0.78 


0.89 


5.3 


(37r) dir 


2.5 


2.6 


0.40 


2.1 


0.023 


0.91 




5.1 


5.3 


0.76 


8.0 








(47r)dir 


12.5 


13.0 


0.13 


3.4 


5.6-10-6 


1.1-10-3 


iru 


0.6 


0.62 


0.75 


0.93 








27rp 


3.6 


3.7 


0.074 


0.55 








2p 


0.9 


0.93 














(57r)dir 


31.0 


31.1 


0.032 


2.1 








27r77 


2.0 


2.1 


0.16 


0.66 








2TrL0 


9.2 


9.5 


0.070 


1.33 










2.3 


2.1 














(67r)dir 


17.2 


17.8 


4.2-10-'^ 


0.15 








UJTj 


1.5 


1.6 














2u 


3.0 


3.1 














(77r)dir 


5.9 


6.1 


1.9-10-3 


0.023 










1.0 


1.0 















are strongly suppressed or even completely closed due to reduced Q values. For instance, 
the (57r)dir channel, most important in the vacuum, is suppressed by a factor ~ 15 in the 
case SI. All channels with heavy mesons r],p,uj are closed in the case SII. 

Table VIII presents the A'"A'" — > nn widths for different total number of pions in the 
final state n . These widths are calculated by summing partial widths of channels with 
the same pion multiphcities in Table VII. The total annihilation width is obtained by 
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TABLE VIII: The widths of NN —>■ mr channels as well as the total width of N annihilation 
(all in MeV) for different parameter sets introduced in Eqs. (71)-(73). 



n 


SO 


SI 


SII 


2 


0.4 


0.8 


5.3 


3 


7.9 


10.1 


0.9 


4 


18.2 


4.8 


10-3 


5 


46.1 


4.0 





6 


22.5 


0.2 





> 7 


8.9« 


0.02* 





total 


104 


19.9 


6.2 



"Obtained as difference between J2^0c — 104 MeV (see Eq. (70)) and the total width of NN nn 

c 

channels with n < 6 . 

Only n = 7 contribution included. 



summing all partial contributions, 

r^E^c- (74) 

c 

Based on results presented in Table VIII, we conclude that the N annihilation width in 
the medium can be suppressed by large factors ~ 5 (SI) or even 15 (SII) as compared to 
the naive estimate of the case SO. Corresponding life times are 

T = ^ ~ 1.9 (SO), 9.9 (SI) , 32 (SII) fm/c , (75) 

Life times in the range 10-30 fm/c open the possibility for experimental studies of bound 
antinucleon-nucleus systems. 

Even larger suppression of annihilation may be expected for nuclei containing anti- 
hyperons. Unfortunately, due to absence of detailed data on BN annihilation, it is not 
possible to perform analogous studies for S = A, E , . . . Available data [54, 55] on Ap 
annihilation show that a-^'^ ^ cr^f"^ at piab ~ lOGeV/c. On the other hand, due to the 
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100 




FIG. 30: Probabilities of the NN mr annihilation as a function of the pion multiplicity n. 
Dotted lines corresponds to the Gaussian fit [39] of experimental pp data. 

appearance of a heavy kaon already in the lowest-mass annihilation channel, AA^" — > ttK, 
one may expect significant in-medium suppression of this channel, contrary to the reaction 
NN — > 27r, which is enhanced as compared to the vacuum (see Table VIII). 

It is worth to mention another qualitative effect, which accompanies annihilation of 
antibaryons in nuclei: the strong modification of the distribution P„ in the number of 
secondary pions, n , as compared to the BN annihilation in vacuum. These distributions 
are calculated as P„ = ^bn-^ti-k/^ presented in Fig. 30 for B — N. The dashed curve 
corresponds to the SO data from Table VIII [79]. In this case the maximum is at n = 5 
and the shape is not far from Gaussian (dotted line) . It is clearly seen that the maximum 
of Pn shifts to smaller values, n = 3 and 2, in matter with densities po (SI) and 2po (SII). 
This effect may be used to select events with formation of superbound nuclei containing 
antinucleons, e.g. by rejecting events with more than two soft pions in the final state. 
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E. Annihilation of antibaryons in finite nuclei 



In the end of this section we estimate the hfe times of bound antibaryon-nucleus 
systems with respect to annihilation. Since the antibaryon is locahzed in a central core 
of the nucleus, our former assumption of a homogeneous antibaryon density is not valid. 
The minimal energy available for annihilation is now given by 

^/sl = Q = mm{E^ + En) = m-^ + tun - - Em , (76) 

where B^ and B^ are the binding energies of the annihilating partners. The annihilation 

rates can be calculated by averaging the local expression (56) over the volume where 
the antibaryon wave function ip^{x) is essentially nonzero. Taking into account that the 
antibaryon density distribution is normalized to unity, after a simple calculation we 
get 

r = J2^c, r, = roScAc(s*)JBiv (77) 

c 

Here we have introduced the overlap integral appropriate for localized density distribu- 
tions, 

J r^yPsBPSN^ (78) 

We have performed numerical calculations for the system using the NLZ2 and 
NL3 parameter sets with two choices of antiproton coupling constants corresponding to 
^ = 1 and ^ = 0.5 . The results are summarized in Table IX. One can notice strong 
model dependence of JpN which is mainly caused by its high sensitivity to the in-medium 
effective masses [80]. As one can see from Fig. 9, the (anti)nucleon effective mass in the ^^O 
system becomes rather small in the central region, especially for the NL3 parameter set. 
This explains large value of JpN in this case. 

The appearance of effective masses in the overlap integral (78) is an artifact of quasi- 
classical approximation used in our study of the in-medium annihilation. For its validity 
the effective potential should vary slowly over a Compton wave length of corresponding 
particles, A = 1/m* (j = B, N) . The approximation breaks down when m* become small. 
Therefore, one can use Eq. (78) only as a rough estimate [81]. 
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TABLE IX: Characteristics of antiproton annihilation in the ^^O system calculated within the 
NLZ2 and NL3 models.The third and fifth columns shows the results of calculation with reduced 
antiproton couplings = 0.5). 





NLZ2 


XL3 




1 


0.5 


1 


0.5 


Bn (MeV) 


88 


83 


113 


102 


Bp (MeV) 


897 


500 


1079 


583 


(^iViev ) 


892 


1294 


685 


1192 




0.96 


0.99 


0.92 


0.98 




0.13 


0.40 


0.046 


0.32 




8.4-10-3 


0.13 


2.9 • 10-"^ 


0.078 


c 


0.81 


9.4 


0.47 


6.9 


JpN 


10.1 


8.5 


46.3 


33.1 


r (McV) 


8.5 


83 


23 


236 


T (fm/c) 


23 


2.4 


8.7 


0.84 



" Calculated as average between single-particle binding energies Bp and Bn (see Fig. 10). 



The estimated total annihilation widths for the ^^O system are also presented in Ta- 
ble IX. In the case of pure G-parity (^ = 1) they arc in the range 9-23 MeV, depending 
on the RMF parametrization. The corresponding life times are 9 - 23 fm/c, in agreement 
with the values presented in Eq. (75). One can indeed talk about the delayed annihila- 
tion due to in-medium effects. These results are very promising from the viewpoint of 
experimental observation of deeply bound antibaryon-nuclear systems. 

However, much larger annihilation widths are predicted by calculations with reduced 
BN couplings (^ = 0.5). This is a consequence of increased Q values and larger number 
of open annihilation channels as compared to the case ^ = 1 . The corresponding life 
times of about 1-2 fm/c are perhaps too short for pronounced observable effects. We 
present here these results to demonstrate uncertainties in our present knowledge of the 
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in-medium annihilation. 

F. Multi— nucleon annihilation 

In addition to the in-medium effects considered above, there are several other processes 
which can inffuence the antibaryon annihilation inside the nucleus. Basically, one should 
consider two types of processes. Processes of the first type are of the long-range nature. 
They include rescattering and absorption of primary mesons on the intranuclear nucleons. 
Since primary pions produced in the annihilation have energies in the A-resonance region, 
they can effectively be absorbed in the chain of reactions: Tr + iV-^'A, A + A^^A^ + iV. 
As a result of meson-nucleon and nucleon-nucleon rescatterings, multiple particle-hole 
excitations will be produced. The nucleus will eventually heat up and emit a few nucleons. 
These processes can be well described by the intranuclear cascade model (see Ref. [56]). 
We expect that they cannot change noticeably the life time of the trapped antibaryon. 

Processes of the second type may, in principle, significantly affect the annihilation 
probability. They include new annihilation channels which arc not possible in vacuum, 
for instance, the emission and absorption of a virtual meson. Since a virtual particle can 
propagate only within a Compton wave length (1/m^ ~ 1.4 fm for pions), this process 
should be of the short-range nature. In other words, the recoiling nucleon should be 
very close to the annihilation zone which is characterized by the radius of about 1 fm [3] . 
Thus, in addition to the annihilation on a single nucleon, considered so far, there exist new 
annihilation channels involving two and more nucleons. It is customary to classify these 
annihilation channels by the net baryon number of a combined system, i.e. B = 0, 1, 2,... 
Obviously, the channels with B ^ must contain not only mesons but also baryons in the 
final state. Probably, the most famous process of this kind was proposed by Pontecorvo 
many years ago [57]: 

p + d^p + 7r~ . (79) 

It has only one pion and one nucleon in the final state {B ~ 1). Despite of a seemingly 
simple final state (two charged particles) , up to now this process has not been convincingly 
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identified experimentally. One may expect that an analogous process with B — 1, 

B + (NN) ^ M + N, (80) 

can also occur on a correlated 2-nucleon pair in heavier nuclei. The next most interesting 
reaction of this kind would be the annihilation on a 3-nucleon fluctuation: 

B + {NNN)^N + N. (81) 

The relative probability of multi-nucleon annihilation can be estimated on the basis of 
a simple geometrical consideration. Let us assume that the annihilation proceeds through 
an intermediate stage when an "annihilation fireball" [58] is formed. It is widely accepted 
(see e.g. Ref. [3]) that within this fireball baryons and antibaryons are dissolved into 
their quark-antiquark-gluon constituents. The radius of this fireball, i?ann, can be found 
from the cross section of the inverse reaction p + p^p + p + N + N. This cross section 
is certainly only a fraction of the inelastic pp cross section, crf^, which is about 30 mb 
in the vicinity of the pp production threshold. Assuming that af^ ~ 27rRlj^^ we get 
-Rann — 0.8 fm, which is in agreement with other estimates [3, 59]. 

Now we can estimate the average number of nucleons, which are present in the anni- 
hilation zone around an antibaryon 



n = PN{r)-—Rt^,. (82) 



Taking pNij") — 2po= 0.3 fm~^, one obtains n ~ 0.6. Finally, we assume that the actual 
number of nucleons present in the annihilation fireball is distributed according to the 
Poisson law, P„ — n"exp [—n)/n\. This gives the probabilities of different channels, 

Po = 0.55, Pi = 0.33, P2 = 0.10, P3 = 0.02, P4 = 0.004, ... (83) 

The relative probability of multi-nucleon annihilation channels is now given by 

£^^0.10+0.02+. 

Pi 0.33 ^ ' 

Therefore, we expect that the channels with P > may lead to 40% reduction of the 
antibaryon life times in nuclei, as compared to estimates including only the P = channel. 
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Experimental information on multi-nucleon annihilation channels is very scarce. We 
mention here a recent paper by the OBELIX collaboration [60], where the p annihilation at 
rest in "^He was studied. They have analyzed the annihilation channels with 27r~^ and 27r~, 
with and without a fast proton (piab > 300 MeV) in the final state. The events where such 
a fast proton was present were associated with the annihilation on a 2-nucleon fluctuation 
i.e. with the B — 1 channel. The corresponding branching ratios are 

B{p +^ He ^ 27r+27r-) = (1.42 ± 0.19)% , (85) 
S(p+^ He ^ 27r+27r-pfast) = (0.098 ± 0.02)% . (86) 

Thus, these data indicate that the relative contribution of the B — 1 channel is less 
than 10%. Certainly, more exclusive data are needed to assess the importance of multi- 
nucleon annihilation channels. These processes, in particular the Pontecorvo-like reac- 
tions, are interesting by themselves, irrespective of their contribution to the total anni- 
hilation rate. They may bring valuable information about the physics of annihilation in 
nuclei. 

VI. FORMATION IN pA REACTIONS 

In this section we present estimates of formation probability of deeply bound 
antibaryon-nuclear systems in pA reactions. It is well known from pp experiments that 
the annihilation cross section o'ann(\/s) very large at low energies [48] (see its low-energy 
parametrization given by Eq. (46)). By this reason slow antiprotons in pA interactions are 
absorbed at far periphery of nuclear density distribution. This situation can be avoided 
by using high-energy antiprotons, whose annihilation cross section is strongly reduced so 
that they can penetrate deeply into the nucleus. 

As most promising from the experimental point of view we consider the following 
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reactions 



p + N ^N,,,^ + N + 7T, £;th = 0.3GeV, (87) 

p+iV^ Asiow + A + TT, Eth = 1.13GeV, (88) 

p+N ^Ki^^ + N + K, £;th = 1.58 GeV, (89) 

p + N^N,,,^ + N + N,,,^ + N, £;th = 5.63GeV, (90) 

p + N ^N,i,^ + N + K,,^ + A, Eth = 7.11GeV, (91) 



where N stands for nucleons (p, n) , E'th is the threshold energy for a corresponding chan- 
nel. Here slow antibaryons have a chance to be trapped by a target nucleus leading to 
the formation of a deeply bound antibaryon-nuclear system. The fast reaction products 
{N, TT, A, K) can be used as triggers. 

One may wonder, why do we need antiproton beams? Wouldn't it be easier to use for 
the same purpose much cheaper and widely available proton beams? The reason is that 
the threshold energy for the baryon-antibaryon pair production in pp collisions is very 
high, at least 5.6 GeV for producing a pp pair. This means that the reaction products are 
fast in the lab frame where the target nucleus is at rest. In this situation the capture of 
an antibaryon by the target nucleus is strongly suppressed. In the case of the antiproton 
beam the corresponding threshold energies are quite low, as in reaction listed above, so 
that the antibaryon capture becomes in principle possible. It is interesting to note that 
antiproton beams open unique possibility to produce simultaneously two antibaryons (the 
last two reactions above) which can then form a bound state like an antideuteron. 

In general the formation probability of a deeply-bound antibaryon-nuclear state can 
be expressed as 

-Pform ^cent ' ^surv ' ^stop ' ^capt ) (^^) 

where Wcent is a fraction of central events selected in a given experiment (typically 
Wccnt — 10%), "if surv IS the autiprotou survival probability. The last two factors in Eq. (92) 
give, respectively, the stopping and capture probabilities. For estimates we assume that 
an antiproton penetrates to the center of a nucleus, i.e. traverses without annihilation a 
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distance of about R (the nuclear radius) in target matter of normal density, 

Wsurv = exp (-po(7ann-R) ■ (93) 

To be captured into a deep bound state an incident antiproton must change its energy 
and momentum in, at least, one inelastic collision inside the nucleus. This can be achieved, 
for instance, by producing a fast pion or kaon carrying away excessive energy and mo- 
mentum. The probability of such an event is denoted by Wstop in Eq- (92). Obviously, 
the energy loss should be equal to the energy difference of the initial and final nuclei. The 
final B momentum should be comparable with the momentum spread of its bound-state 
wave function. One can estimate this momentum spread as Ap ~ ir/Rp, where Rp is the 
rms radius of the antibaryon density distribution. Prom numerical calculations presented 
in Sect. IV we find ~ 1.5 fm so that Ap ~ 0.4 GeV. 

The probability of a single inelastic pN collision can be estimated by assuming Pois- 
sonian distribution in the number of such collisions, Wn = n"" exp (— n)/n!, so that 
wi — nexp (— n). The mean number of inelastic collisions on a distance r = R is 

R 

n = — = poo-in^ , (94) 

where Ain is the mean free path between inelastic collisions and is the total pN inelastic 

cross section. 

In fact only a small fraction of inelastic collision leads to a desired energy-momentum 

loss. This fraction can be calculated from the differential inelastic cross section, for 
_ d'^p 

the reaction pN BX. Explicitly we define Wstop — ■ wi^ss with 

wio,s = w{p^^ <Ap)^— / dapj^_^BX ■ (95) 

Experimental data on the reaction pN BX in a GeV energy region are quite scarce. We 
have found only one paper [61] where detailed results from the bubble chamber experiment 
at CERN were reported. The differential cross sections for the reaction pp AX were 
measured for several p momenta from 3.6 to 10 GeV/c. Using these data we estimate 
w{p^^ < Ap) in the range 10~^ -j- 10~^ for p incident momenta 3.6-=-12 GeV/c. The 
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integrated cross section for the reaction pp — > AX at 3.6 GeV/c is approximately 0.5 mb 
as compared to about 50 mb for the total inelastic pp cross section. The ratio of these 
two values, ~ 10"^, gives the probability of the p — >■ A conversion, which is small in 
accordance with the Zweig rule. In considered region of bombarding energies we can 
obtain 'w{pp^ < Ap) simply rescaling w{p^^ < Ap) by a factor ~ 10^. 

TABLE X: Formation probabilities and reaction rates for producing antibaryon-nuclear systems 
in several typical reactions. The values of Pform are obtained by multiplying probabilities from 
the first three columns with Wcentuicapt , where Wcent = if capt = 0.1 . The reaction rates are 
obtained as y = Pform ■ L ■ a^eom with L = 2- 10^^ cm-^s"^, a^eom = 0.27 and 1.1 b for ^''O and 
^^^Bi, respectively. 



piab (GeV/c) 


reaction 




Wl 


W^loss 


Pform 


Y (s-i) 


0.8 




0.03 


O.IO** 


0.1 


3 • 10-6 


160 


3.6 




0.26 


0.36 


10-2 


9 • 10-6 


510 


3.6 




0.26 


0.36 


10-^ 


9 • 10-s 


5 


3.6 


|5+209Bi^208pb + „ + ^+ 


0.07 


0.29 


10-2 


2 • 10-6 


440 


3.6 


p+209Bi^208pb+p + ;^0 


0.07 


0.29 


10-^ 


2 • 10-^ 


4 


10 




0.53 


0.37 


10-3 


2 • 10-6 


110 


10 




0.53 


0.37 


10-5 


2 • 10-^ 


1 


10 


p+209Bi^208pb + „ + ^+ 


0.22 


0.21 


10-3 


5 • lO-'^ 


110 


10 


p+209Bi^208pb+p + ;^0 


0.22 


0.21 


10-5 


5 • 10-9 


1 



° In this estimate we use the value a\n 3 mb motivated by measurements for the reactions p N — > tt 
at Ef'' = 300 MeV reported in Ref. [62]. 

The last factor in Eq. (92), Wcapt, is the most uncertain quantity. It is determined by 
the matrix element between the initial, plane wave state and the final, localized B state. 
Moreover, because of the nuclear polarization effects, in particular, the reduced effective 
mass and finite width, the bound antibaryon is off the vacuum mass shell. The realistic 
calculation of tfcapt will require a special effort which is out of the scope of the present 
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paper. For our estimates below we take a fixed value tWcapt = 0.1. 

In Table X we present resulting formation probabilities for several typical reactions. 
One can see that Pform is on the level of 10"''-=- 10~^ for bound p + A systems. In this Table 
we also give reaction rates calculated for the parameters of antiproton beams planned at 
the future GSI facility [63] (the beam luminosity is L = 2 • 10^^ cm~^-s~^, antiproton 
energies from 30 MeV to 15 GeV). We see that expected reaction rates are in the range 
from tens to hundreds desired events per second. Such rates seem to be well within the 
present detection possibihties. 

Let us consider now two examples which clarify the reaction kinematics. For simplicity 
we use plane-wave functions labeled by the particle momentum p. We are interested 
in a reaction which leads to the final B trapped in a deeply bound state with energy 
Ej — ^ {rrij — Sj)"^ + Pj + Vj, where j — N,A, . . . For instance, consider the reaction (87) 
where an antiproton with initial energy E'beam is captured in a target nucleus (^4, Z) after 
colliding with a bound nucleon N with energy E]\r. If the recoiling nucleon, say a neutron, 
leaves the nucleus the energy of the final nucleus p{A — 1, Z) is changed by 

A£; = 4-K- [^{rrip-Spy+p^' + Vp^ - (mn - + pI + V^) , (96) 

where Pn is the initial neutron momentum and is the final antiproton momentum, both 
defined inside the nucleus. Assuming that p^ ~ p„ and ~ S'„, Vp — —Vp (G-parity) we 
find that the binding energy of the final nucleus is changed by AB — —AE ~ (Vp + Vn). 
Neglecting small recoil effects, we now estimate the energy of the emitted pion as 
Et^ ~ -E^beam + AB . So, the observation of a pion with energy exceeding the incident 
antiproton energy would provide a strong evidence in favor of the p-nuclear bound state 
formation. Such pions can be well separated from the annihilation pions which are pro- 
duced at a later stage of the reaction. Due to the expected large energy gain such non- 
annihilation pions can be produced even at subthreshold antiproton energies. 

Analogous consideration can be done for the reaction (89) leading to the formation of 
an A- nuclear bound system, j^lA — 1, Z). The energy difference between the initial and 
final nuclei is in this case 

AE ~ (mA - m^) - {Sj, - S^) - {V^ + V^) . (97) 
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With Sa ^ Sn and Va ~ Vn this energy difference is negative. The energy balance 
requires that the energy -Ebcam — AE" is carried away by an emitted antikaon which, 
therefore, should be fast in the lab frame. We propose to use fast pions and antikaons as 
triggers for selecting desired events. 

It is necessary to emphasize that the antibaryon capture reactions discussed above are 
relatively fast as compared to the characteristic time of target nucleus rearrangement, Tnuc- 
Therefore, initially antibaryons will occupy "non-modified" levels which are predicted for 
ordinary nuclei. Rearrangement of nuclear structure, in particular, the creation of a local 
compression, as discussed in Sect. IV, will take a longer time. One can estimate this 
time as 

''"nuc — Rcore/ C-s i (^^) 

where i?core is the radius of a compressed core and is the sound velocity. Taking 
i?corc — 1-5 fm and Cg = 0.2 c we get Tnuc — 7.5 fm/c. In the light of estimates presented in 
Sect. V, this time is of the same order as the predicted life time of antibaryon-doped nuclei. 
Better understanding of the rearrangement process can be achieved only by carrying out 
dynamical simulations on the basis of the TDHF or molecular dynamics approaches. 

Finally we estimate the temperature of an antibaryon-nuclear system which might be 
expected after the rearrangement of its structure. As one can conclude from Fig. 18 after 
the rearrangement the antiproton binding energy increases by about 220 MeV in ^^0 
and 400 MeV in |He. Assuming that this binding energy gain AE is transferred into 
heat, we can estimate the temperature of the final system using the Fermi gas formula 
AE — aT"^ , where a ~ ^4/10 MeV is the nuclear level density parameter. For this gives 
T ~ 12 MeV which is much lower than the temperatures ~ 100 MeV associated with high- 
energy nuclear collisions. Therefore, we can indeed speculate about cold compression of 
nuclei induced by antibaryons as an alternative to the shock-wave compression in heavy- 
ion collisions [64]. In the latter follows from the Hugoniot adiabate, compression 
of nuclear matter is always accompanied by its strong heating. 
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VII. OBSERVABLE SIGNATURES OF BOUND ANTIBARYON-NUCLEAR 

SYSTEMS 



If bound antibaryon-nuclear systems exist and live long enough they can manifest 
themselves in several ways. In this section we discuss their possible signatures. 

A. Transitions from atomic to nuclear states 

If a slow antiproton is first captured on a Coulomb orbit forming an antiproton atom, 
it can later on make transition on a deep nuclear bound state. This process was first 
proposed in Ref. [5] where antiproton-nuclear bound states were studied within the non- 
relativistic optical model. Such a transition will be accompanied by the emission of a 
monoenergetic photon, pion or kaon, depending on the antibaryon type. According to 
our calculations, many discrete states may exist in the relativistic antiproton potential. 
Their binding energies and thus, the transition energies may vary from a few tens to 
several hundred MeV. The probability of such transitions is determined by the matrix 
element of the corresponding transition operator between the Coulomb and nuclear an- 
tiproton wave functions. As shown in Ref. [5], the partial width associated with such 
transitions could be as small as 10"^ — 10~^ of the total width of the atomic level. The 
total width is obviously determined by the annihilation on the nuclear surface. There- 
fore, a special effort should be made to find this signal in the huge background of direct 
annihilation events. 

B. Super— transitions from upper to lower well 

As pointed out earlier, high-energy antiprotons have a better chance to penetrate deep 
into the nucleus. Under certain conditions discussed earlier they can be trapped into a 
deep bound state. We have demonstrated in Sect. V that the annihilation can be delayed 
in this case due to the reduction of energy released in this process. An antibaryon sitting 
in the nucleus can be viewed as a hole in otherwise filled Dirac sea. This state represents 
a strong excitation of the nucleus which soon or later will decay. Among other decay 
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TT, K,7 




FIG. 31: Annihilation of antibaryon-nuclear system as a "super-transition" between the Fermi 
and Dirac seas. 

modes the most interesting is a super-transition when a nucleon from a discrete level of 
the Fermi sea jumps into the hole state in the Dirac sea. This can be achieved by emitting 
a single photon, pion or kaon with energy of about Q = Ej^ + E-j^ and isotropic angular 
distribution in the nucleus rest frame. Such process is forbidden in vacuum by the energy- 
momentum conservation laws. But in the considered case the recoil momentum can be 
carried away by the residual nucleus. This is analogous to the Pontecorvo-like reactions 
discussed in Sect. V. This process is illustrated in Fig. 31. The appearance of relatively 
narrow lines with energies of about 0.5 — 1 GeV and width AE ~ 20 ^ 50 MeV in the 
spectrum of secondary photons or mesons would be a direct signature of the deeply-bound 
antibaryon-nuclear states. 

C. Explosive multifragmentation 

Another signal may come from the explosive disintegration of a compressed nucleus 
after the antibaryon annihilation. When the driving force for the compression disappears, 
the nucleus will expand as a compressed spring. As a result of the collective expansion the 
nucleus will be torn apart into fragments, as illustrated in Fig. 32. This process can be ob- 
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FIG. 32: Multifragmentation (right) of a nuclear remnant formed after the antibaryon annihi- 
lation in a bound B-nuclear system (left). 

served by measuring collective velocities of these fragments. Multifragmentation of nuclei 
induced by high-energy antiprotons has been studied earlier [65], but with limited statis- 
tics. So far only minimum bias events have been analyzed. In this case the distribution 
of fragments follows closely predictions of the statistical multifragmentation model [56] 
where no collective effects are included. We emphasize again that a dedicated study with 
proper triggering is needed to find explosive events. Modern experimental technics allow 
to unambiguously distinguish between purely thermal and flow-driven multifragmenta- 
tion. 

D. Multi— quark— ant iquark clusters 

It is interesting to look at the antibaryon-nuclear systems from a somewhat different 
point of view. An antibaryon implanted into a nucleus acts as a strong attracting center 
for nearby nucleons. Due to uncompensated attractive force these nucleons acquire accel- 
eration towards the center. As the result of this inward collective motion the nucleons will 
pile up and produce a compression zone around the antibaryon. If such a process were 
completely elastic, it would look like a monopole-type oscillations around the equilibrium 
configuration found by solving static equations (see Sect. IV). In this dynamical process 
even stronger compression can be reached as compared with the one predicted by this 
static configuration. The maximum compression is achieved when the initial potential 
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FIG. 33: Schematic picture of deconfined quark-antiquark droplet formed in interaction of an 
antiproton with the ^He nucleus. 

energy generated by the antibaryon is transformed into the compressional energy. Simple 
estimates show that local baryon densities up to 5po may be obtained in this way. It is 
most likely that the deconfinement transition will occur at this stage and a high-density 
cloud containing an antibaryon and a few nuclcons will appear in the form of a multi- 
quark-ant iquark cluster. One may speculate that the whole "^He or even ^^O nucleus 
can be transformed into the quark phase by this mechanism. This process is illustrated 



As shown in Refs. [33, 34] and shortly discussed in Sect. Ill (see Fig. 4), an admixture 
of antiquarks to cold quark matter is energetically favorable. The problem of annihilation 
is now transferred to the quark level. But the argument concerning the reduction of 
available phase space due to the entrance-channel nuclear effects should work in this case 
too. Thus one may hope to produce relatively cold droplets of the quark phase by the 
inertial compression of nuclear matter initiated by an antibaryon. A similar mechanism 
for producing large quark bags was proposed in Ref. [58]. There the idea was that the 
annihilation fireball moving through the nucleus can grow by absorbing nucleons on its 
way. 

Thus, the annihilation process in nuclear environment may serve as a breeder for cre- 
ating new multi-quark-antiquark structures. We believe that among the decay modes 
of multi-quark-antiquark systems, like the one depicted in Fig. 33, not only usual 
hadrons, but also exotic states such as pentaquarks {uuddd, uudds, . . .) or heptaquarks 
{uuuuddd, uuuddd, . . .) can be present. Moreover, more exotic states like baryonium 



in Fig. 33. 
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(3g3g), deutronium {6q3q), tritonium (9^3^) and even helionium {12q3iq , see Fig. 33) 
can also be searched for. Discovery of such states would significantly extend our knowl- 
edge about the quark-gluon structure of matter. Such a study will be complimentary 
to the active search for exotic multi-quark systems stimulated by the discovery of the 
pentaquark [66]. 

VIII. SUMMARY AND DISCUSSION 

In this paper we have demonstrated that antibaryons are of significant interest for 
nuclear physics. They can provide new and important information about the nature of 
strong interactions and structure of the vacuum. As well known, the relativistic nuclear 
models predict deep antibaryon potentials in nuclei generated by coherent action of mean 
scalar and vector fields. To verify this picture we study a real antibaryon bound in the 
nucleus instead of considering filled states of the Dirac sea. On the basis of the RMF 
model we performed detailed calculations of such bound antibaryon-nuclear systems tak- 
ing into account the rearrangement of nuclear structure due to the presence of antibaryons. 
The self-consistent calculations lead to stronger bindings as compared to previous stud- 
ies where the rearrangement effects were ignored. What is even more important, our 
calculations predict strong local compression of nuclei induced by the antibaryon. This 
opens a principal possibility of producing relatively cold and compressed nuclear matter 
in the laboratory. In contrast, in high-energy nuclear collisions the compression is always 
accompanied by a strong heating of nuclear matter. 

Our second goal in this paper was to deceive a common delusion that the antibaryon an- 
nihilation in nuclei is so strong that it will shadow all other processes. We have performed 
detailed calculations assuming that annihilation rates into different exclusive channels are 
proportional to the available phase space. Due to reduced effective masses and (partial) 
cancelation of vector potentials, the energy available for annihilation of slow antibaryons 
(the reaction Q-value) is significantly reduced as compared to the minimum vacuum value, 
rriN + iTi-g. Then many channels of in-medium annihilations are simply closed. This may 
lead to a dramatic reduction, by factors 5-^20, of the total annihilation rate. We estimate 
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life times of deeply-bound antibaryon- nuclear systems on the level of 2-=-20 fm/c that 
makes their observation feasible. This large margin in the life times is mainly caused by 
uncertainties in antibaryon coupling constants as well as the overlap integrals between 
the antibaryon and nucleon scalar densities. We have also analyzed multi-nucleon anni- 
hilation channels (Pontecorvo-like reactions) and found their contribution to be less than 
40% of the single-nucleon annihilation. 

We believe that bound antibaryon- nuclear systems can be produced by using antiproton 
beams of multi-GeV energy, e.g. at the future GSI facihty. Since the annihilation cross 
section drops significantly with energy, a high-energy antiproton can penetrate deeper 
into the nuclear interior. Then it can be stopped there in an inelastic collision with a 
nucleon, e.g. via the reaction A{p, Ntt)pA', leading to the formation of ap-doped nucleus. 
Reactions like A {p, N K)-^A' a,ndA{p, A7r)j;A' can be used to produce a A-doped nucleus. 
Fast mesons, nucleons or lambdas can be used for triggering such events. Our estimates 
of the formation probability in a central pA collision give the values 10~^ — 10~^. With 
the antiproton beam luminosity of 10'^^ cm~^ s~^ planned at GSI this will correspond to 
the reaction rate from tens to hundreds of desired events per second. 

We have proposed several observable signatures which can be used for detection of 
antibaryon-nuclear bound states. They include exotic annihilation channels with emission 
of a single photon, pion or kaon; explosive multifragmentation of nuclei, and formation 
of multi-quark-antiquark clusters. The possibility of using antibaryon annihilation in 
nuclei to produce droplets of relatively cold quark matter is interesting by its own. This 
mechanism will work irrespective to the existence of antibaryon-nuclear bound states. 

There remain many problems to be studied. The list of most urgent questions includes: 

• how G-parity works in a many-body system and what are the antibaryon couplings 
to the meson fields? 

• what is the role of exchange terms and dispersive corrections? 

• how large is the unphysical self-interaction of antibaryons within the RMF approach? 

• how to extend the RMF approach to include self-consistently imaginary contribution 
to the antibaryon self energy? 

• how accurate are quasiclassical estimates of life times of bound antibaryon-nuclear 
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systems? 

• how these hfe times are affected by in-medium modifications of secondary mesons? 

In addition, the nuclear rearrangement dynamics after the antibaryon capture should be 
studied on the basis of a dynamical approach hke the TDHF model. Transport calculations 
of pA reactions are needed to check our simple estimates of the formation probabihties 
and annihilation rates. We are planning to address these questions in the future work. 

Our most general conclusion is that the antibaryon-nuclear physics is a broad field of 
research which is not yet explored sufficiently well. It may bring about new interesting 
phenomena ranging from unusual annihilation channels to exotic antibaryon-nuclear or 
multi-quark-antiquark bound states. This is a nonperturbative domain of strong interac- 
tions and therefore the progress in this field can be achieved only by close cooperation of 
theory and experiment. 
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to the local density i.e. requires an additional nucleon in the annihilation zone. Moreover, 
the s-wave contribution for pions and p-mcsons is proportional to the isovector density 
which vanishes in the isosymmetric matter. We classify such processes as multi-nucleon 
annihilation and consider them in Sect. VF. The long-range contributions are generated by 
the nucleon-hole and resonance-hole diagrams which are characterized by the energy scale 
of about 100 MeV. Therefore, they can only affect propagation of mesons at large distances 
(exceeding about 2 fm) or at long times. This is why we believe that these in-medium 
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[77] Here and below we consider the static B + N matter in its rest frame. 
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should be compared with the rearrangement time. 

[79] To obtain P„ at pjv — , one can use the SO results from this Table since in this case pN 
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•^nnP^/ ~ ■'■^•^ obtained from the previous estimates for infinite matter in the case SII. 

[81] In fact, the overlap integral is proportional to / = (^'^~^^ where Ej = ^Jra*^ ^ (r) + is 
the kinetic energy of the j-th particle and angular brackets denote averaging over single- 
particle wave functions of nucleons and an antibaryon. Due to a finite momentum spread 
of the wave functions / is always finite even if m*{r) vanishes. This is why we think that 
the quasiclassical expression (78) which contains m* in the denominator overestimates the 
overlap integral and thus, the annihilation width of a bound antibaryon. 
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